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ABSTRACT 


The spectral distribution and spectral density of 
electromagnetic modes are perceived to be of value in the 
design and analysis of rectangular multimode cavities. 
However, these mode distributions have never been determined 
correctly in the engineering literature. Serious errors 
have been discovered in both numerical and analytical mode 
distribution calculations which have survived, apparently 
undetected, partly due Ee a lack Of experimental 
measurements. The nature of the errors suggests that 
misunderstandings exist regarding the partial distributions 


of TM and TE modes. 


Agsbrief introductions isSeugiven Jewhich corrects a 
deficiency in the engineering literature regarding 
references to fundamental theoretical WORke The 
introduction is followed by a review of the essential 
physics of electromagnetic boundary value problems in 


rectangular geometry. 


A numerical algorithm is presented which enumerates TM 
and TE modes’ separately, and thereby yields the correct 
number of modes in an empty rectangular cavity for any 
finite bandwidth. This algorithm is used to generate 
corrected versions of previously published cavity design 


tables which were based on faulty algorithms. Numerical 
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mode distributions are computed for a number of hypothetical 
cavities which illustrate some important characteristics of 
these distributions not previously considered. These 
characteristics affect the utility of mode distribution 


calculations in the design of microwave heating cavities. 


The first experimental measurements of mode 
distributions are reported for an empty rectangular 
laboratory cavity which are in good agreement with computed 
distributions obtained using the correct algorithm. 
Additional results for the same cavity fitted with a mode 
stirer indicate that the mode distributions are 


Substantially perturbed by such loading. 


Certain versions of the scalar asymptotic mode 
distribution formulae, which are valid for acoustic modes, 
have been incorrectly applied to microwave cavity problems 
in engineering. The correct formulae are derived by 
considering the number of lattice points in an ellipsoid in 
p-space. Approximate eigenfrequency distributions obtained 
from these formulae are compared with the numerically and 
experimentally determined distributions. The results of 
these comparisons verify the analytical formulae at 
microwave frequencies and, therefore, are of interest to 
both physicists and engineers. The accuracy of these 
formulae are discussed in terms of their role in the design 


and analysis of microwave power applicators. 
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Cavity. 

average (signed) error and average absolute error 
respectively, between exact and approximate mode 


distrrbutvons 


CHAPTER 1. INTRODUCTION 


The first section of this chapter is intended to 
provide a qualitative description of the nature of 
elgenvalue spectra, and to introduce functional definitions 
of the two types of mode distribution. Proper mathematical 
definitions will be given in Chapter 2. The next two 
sections comprise a brief overview of some of the important 
existing literature, which iS necesSary in order to 
understand the origin of the problems addressed in later 
chapters. The final two sections describe the objectives of 


the present work and define the scope of the thesis. 
1.1 Properties of Eigenvalue Spectra 


The electromagnetic forms of the wave equation follow 
directly from Maxwell's equations and are therefore 
inherently involved in virtually every problem in 
electrodynamics. The classical solution to these equations 
in a domain completely bounded by perfectly conducting 
surfaces yields an infinite set of eigenvalues representing 
the allowable modes of oscillation. These eigenvalues are 
quasi-randomly spaced throughout the spectrum and exhibit 
finite degeneracy at each of the eigenfrequencies. For any 
of the simple geometries in ane the wave equation is 
separable, any arbitrary number of the eigenvalues can be 


calculated analytically. However, ash is virtually 
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impossible to predict, by strictly analytical means, the 
exact number of modes that will occur within a given 
bandwidth. Yet this spectral density is of considerable 


importance in many scientific and engineering applications. 


The spectral density is derived from the more 
fundamental spectral distribution of eigenmodes. Simply 
stated, the spectral distribution of the eigenvalues ffor a 
particular bounded domain (e.g. waveguide or cavity) is the 
cumulative total mumber of eigenvalues, or modes, with 
elgenfrequencies not exceeding some upper limit, as a 
pEtinctvon” "eotwethenelimat Eirequencyr Ti the spectral 
distribution "function" is evaluated at two different marker 
frequencies, the difference yields the number of modes in 
the bandwidth defined by these two frequencies. This result 
is equivalent to evaluating a spectral density 
"function"—which describes the number of modes per interval 
of bandwidth—at the same marker frequencies. If the two 
frequencies are allowed to be arbitrary, but their 
separation is held constant, then this differential number 
of modes can be (and usually is) also presented as a 


distribution in the statistical sense. 


In order to simplify the discussions to follow, the 
fermmidensitySuewplliagbe Used pthrotighouteethis thesrsigin 
reference to both the evaluated spectral density (i.e. the 


differential number of modes) and the unevaluated or true 
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spectral density of modest. It should be clear from the 


context which form, if a particular one, is required. 


In many contexts, the discussion applies equally well 
to either the spectral distribution or the spectral density 
of modes. It is therefore convenient to introduce an 
ambiguous term "eigenmode distribution" or simply "mode 
distribution", which could mean either or both of these 


distributions. 


Of course, given the nature of the eigenvalues 
described above, neither the spectral distribution nor the 
Spectral density of modes can, strictly speaking, be true 
functions. The reason for this phraseology will become 
clear when the exact mode distributions are defined 
mathematically in the next chapter. Also it will be helpful 
to dispense with the semantic distinction when the exact 
mode distributions are compared graphically to the 
asymptotic approximations, which are true functions of 


frequency. 


1.2 Overview of Physics Literature 


There are two methods of determining the exact 


eigenmode distributions for a given cavity: numerical 


+For most engineering applications, the bandwidth must be 
specified in order to extract useful quantitative 
information. 
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computation and experimental measurement. It is important 
to realize that both are relatively modern methodst. But 
near the turn of the century, long before the advent of 
programmable computing machinery and electronic measuring 
equipment, there was a scientific need for an understanding 
of the behavior of eigenvalue spectra [1-4]. Calculus and 
analytical methods were the most powerful theoretical tools 
of the day, and theorists therefore approached the problem 
by studying the asymptotic behavior of eigenvalue spectra in 
the limit of large wave numbers. That is, the average 
cumulative number of modes, as a function of frequency, in 


the high frequency limit. 


The first detailed investigation of this’ kind was 
undertaken more than seven decades ago by Weyl [5-8] who 
considered the mathematical problemt of calculating the 
asymptotic distribution of the eigenvalues of Laplacian 
operators. Weyl obtained the first term of an asymptotic 
expansion of the mode distribution and proved its 
independence of the shape of the domain. Later 
investigators, notably Courant [9], . Courant and. Hilbert 


fl Ole Cardemanaleiin, Pleijel b124431 _¢Brownell [949 15i,8and 


tALthoughsas student » of sHaA. Lorentz, Johanna Reudler, 
actually calculated (by hand) portions of the exact mode 
distributions for a few simple geometries [Leiden 
usSemta wLOneanio | 2h. 

tWeyl at Pst considered some mnonphysical boundary 
conditions which resembled those of the electromagnetic 
Drobicniembutesthis Gidea notasdetractarromethesmathematical 
importance of the work. 
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Agmon [16,17], extended Weyl's results and improved the 
accuracy of the asymptotic expansion. Most of these authors 
considered the special case of electromagnetic boundary 
conditions, hut the emphasis of this early research was more 
on mathematical rigor than physical applications of the 
theory. The references given above are representative but 
by no means constitute a complete list. An excellent and 
thorough review of the important (mathematical) developments 
up to the time of its publication in 1967 is given by Clark 
ite pe 


Some early work which focused more on the 
electromagnetic problem was done by Muller [19], Muller and 
Niemeyer [20], Niemeyer [21,22] and Pathria [23]. It was 
NOteeulbLL wethe early seventies, however, that the 
electromagnetic problem was fully understood theoretically. 
Also by this time, numerical and experimental methods were 
well-developed, which permitted the theoretical work to be 
tested. Some of the important results from this era will be 
discussed in some detail in Chapter 5, but Case and Chiu 
[24:), @BaliansandeBlochs(25)fe@tBalvan Bande@eDuplanticraa[26]F 
Baltes and Kneubuhl [27], Baltes and Hilf [28], Baltes 
[29,30], and Steinle and Baltes [31,32] were eminent 
GOnNtCLIDUtOLS: The theoretical results obtained for the 
electromagnetic case have been applied mainly to infrared 
andmmopticalm (30-37) @¥problemsS@ in@ physics]. Comprehensive 


literature reviews which include both theoretical and 
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applied research on the electromagnetic (and other) problems 
are given in several articles by Baltes, et al. 
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The physics literature cited above is all "well 
connected" in the sense that subsequent papers built on the 
results of previous research and adequate references were 
provided to assist later researchers. However, two papers 
by Bolt [39] and Maa [40], which appeared simultaneously in 
1939, were not quite so "well connected". Bolt and Maa 
treated the scalar (acoustic) problem and derived asymptotic 
formulae for the mode distributions that were presented as 
extensions to those given in the texts by Courant and 
Hilbert [10] and Morse [41]. No other references to prior 
research were given. This apparently singular breech of 
continuity may have contributed to some of the problems 


which later developed in the engineering literature. 


1.3 Overview of Engineering Literature 


The interests of (electrical) engineers are mainly in 
the electromagnetic applications of mode distributions. The 
most common of which being the investigation of the 
characteristics of rectangular cavities used in microwave 
heating systems. Of course, the precise shape of the 


(exact) mode distributions in a loaded cavity are 
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necessarily different from those in an empty cavityf. Hence 
it would seem that, from an engineering point of view, the 
empty-cavity (or intrinsic) mode distributions may or may 
not be important characteristics. But the concensus among 
microwave oven designers appears to be that the intrinsic 
mode distributions are indeed helpful in predicting the 
performance of the loaded cavity [42-49]. However, it is 
interesting to note that these mode distributions have yet 
to be determined correctly—by numerical, experimental or 


analytical means—ain the engineering literature. 


Tet: was not until the early sixties that mode 
distributions were really considered in the context of 
microwave applicator design. Philips Research Labs 
published a set of tables [42] giving the (incorrect) number 
of modes, within the ISM band at 2450 MHz, for several 
rectangular cavities. These tables were later reprinted in 
a book by Puschner [43], accompanied by a discussion of 
exact and approximate mode distributions. Unfortunately, 
the asymptotic formulae quoted by Puschner are also 
incorrect. They are actually the acoustic formulae given in 
the second edition of the text by Morse [41], which 
incorporated the work of Bolt and Maa mentioned previously. 
Puschner simply doubled each of the terms in Maa's formulae 
HInenearlyeallmpractica lmapDl1catlons, mcavityelOcdSeconsise 
of inhomogeneous, irregularly shaped, lossy dielectric 


materials and therefore substantially change the boundary 
value problem. 
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in order to account for the two different wave types. The 
epronslenngebothsethe snumericalandatasi( tables) and the 
analytical formulae will be explained and corrected in later 


chapters. 


James, et al. [44,45] also quoted Maa's formulae in 
articles on multimode theory. These articles have _ been 
cited in several later engineering works, including the 
recent text by Metaxas and Meredith [46] where both [45] and 


[40] are cited. 


Copson and Decareau [47] avoid the mistake of quoting 
Maa's formulae by giving instead only the first term (for 
which the doubling is valid). But this too is misleading 
Since the first term is only accurate for infinitely high 


frequencies. No references are provided. 


The only known attempt (in English) to point out these 
problems is unfortunately buried within an article on 
Microwave applications in China [49], and is itself marred 
bymar ser) ouSs@typographicalserror mnethesmcorrectediespectrad 
distribution formulaft. Again, no references are given which 
direct the reader to the copious physics literature on _ the 
Subject. In fact, ironically, the only direct link between 
the engineering literature and the physics literature is via 
+A more explicit account in Chinese [50] gives the correct 


asymptotic formulae, but does not contain any references to 
the fundamental physics papers. 
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the papers by Bolt and Maa, which have nothing to do with 
the electromagnetic problem, and which do not provide 


adequate references to gain access to the EM work. 
1.4 Objectives 


Partly due to the lack of a theoretical literature 
base, and partly due to the lack of experimental data, 
long-standing errors abound in the engineering literature on 
multimode cavities. These errors have accumulated for over 
two decades without definitive action to correct them. 
Because the errors are of a theoretical nature, 
misunderstandings regarding the physics of mode 
distributions have resulted among engineers working in the 


field of microwave power. 


It iS a major objective of this thesis to draw 
attention to these problems, and to correct the errors by 
presenting the proper theoretical arguments along with the 
pertinent (correct) results. Attention is given to both the 
numerical computation of exact mode distributions, and the 
analytical calculation of approximate mode distributions. 
Also, the bibliographic references bridge a long-standing 
gap between the engineering and physics literature, and 
should provide other scholars and researchers with a 
sufficient literature base to begin more detailed studies of 


this important subject. 
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A second objective is to report the results of 
experimental measurements which substantiate the physical 
theory at microwave frequencies, and which reveal something 
of the nature of mode distributions in the presence of light 
metallic structures (e.g. mode stirrers). These 
measurements were performed on rectangular cavities but the 


technique is applicable to more complicated geometries. 


Amechirdgmob ectivess1Sus to discuss some of the 
characteristics of exact and approximate mode distributions 
which have been discerned from the experimental measurements 
and a number of numerical case studies. These 
characteristics must be understood in order to evaluate the 
roles and the utility of both exact and approximate mode 
distribution calculations in engineering research and design 


applications. 


aoe Scope 


Most of the existing problems in the engineering 
literature involve mode distributions in empty, rectangular, 
ideal microwave cavities. Accordingly, this thesis 
concentrates on electromagnetic modes in perfectly 
conducting rectangular cavities which, for theoretical 
purposes, enclose a free space medium and, for practical 
(e.g. experimental) purposes, are air-filled. Although much 


of the physics of these cavities also applies to other 
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geometries, loading conditions and wave phenomena, these 
cases will not be studied in detail here. Where 
appropriate, generalizations will be mentioned, but not 
discussed unless there are direct ramifications in the field 
of microwave power engineering. Further limitations to the 
scope of certain sections will be given in the text 


pertaining to that material. 
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CHAPTER 2. THEORETICAL BACKGROUND 


It 1s essential to understand the basic physics 
involved in electromagnetic boundary value problems in order 
to correctly obtain and interpret eigenmode distributions. 
It is therefore the purpose of this chapter to review these 
principles and to present, in a mathematical context, some 
important equations used in subsequent chapters. The two 
dimensional eigenvalue problem is introduced in the first 
section which leads to the definition and classification of 
waveguide modes. The analysis is extended in the next 
section to include cavity modes. The eigenfunctions which 
describe the actual field components are presented in _ the 
third section. Rigorous (mathematical) definitions of the 
mode distributions are given in the last section, along with 
a discussion of how and why the partial distributions of TM 


and TE modes are different. 


2.1 Rectangular Waveguide Modes 


Consider a three dimensional semi-infinite domain in 


free space defined, in Cartesian coordinates (x,,X2,X3), by 


Om< mq Ss L, (eal) 
0 < X2 < L2 Gn 2) 
-~e < X3 < +o i233) 


and bounded by perfectly conducting plane surfaces. 
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Electromagnetic fields in this domain must satisfy the 


appropriate form of Maxwell's equations, viz. 


V xX E = -20B (2.4) 
ot 

Vex eH oman D 2D) 
ot 

Vie ceD eas 0 (Zac) 

Vee Bo =a 0 2) 


where B=uoH and D=e,E and, since the domain is enclosed by 


perfect conductors, the appropriate boundary conditions are 


nx E = 0 (2289 
n xX H = J, (2595) 
fie =D" pe (27610) 
n> B= 0 G2 salblp) 


where n is the unit normal vector directed inward (toward 
the free-space medium) from the boundaries. It can easily 
be shown [51-53] that Eqs.(2.4-2.7) are equivalent to the EM 


wave equations 


Vie hee Uo€007E = 0 (52 lez) 
Ot 

V2He= Uo€007H = 0 (2 ees )) 
One 


subject to the same boundary conditions described by 
Eqs.(2.8-2.11). And if the time dependence is assumed to be 
harmonic, then Eqs.(2.12,2.13) further reduce to the vector 


Helmholtz equations 
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V7E + kgE 


i) 
(oo) 


(282) 


V7H + k§éH 


u 
co) 


(255) 


where kg=uo€ow? the infinite-medium, free space wave number. 
It 1S common to view kyo as the modulus of a so-called wave 
vector ko with components in each of the coordinate 


GUrectilons a.Vige 
Ko = k,a, aP k2ae ote k3a3 (23416) 
from which 


[ko]? = kB = k} 


+ 
~ 
INTs) 
+ 
x 
Wn 


625 17) 


and the direction of ko is everywhere normal to the electric 


and magnetic field vectors. 


Due to the infinite extent of the domain “in the +x; 
directions, it is useful (and reasonable) to assume 
propagating wave functions of the form 


-i(wttk3x;3) 


E(x,t) E(x,,x2Je (ar) 


-i(wtt+k3x;3) 


H(x,t) H(x,,xX2)e G2751/9)) 


where the plus/minus Signs indicate forward/backward 
propagating waves. Furthermore, since the boundary 
conditions on E and H are different, the eigenvalues will, 
in general, be different. Thus it is natural to divide’ the 
field solutions “into ttwoledistinct=sclasses: transverse 


magnetic (TM or E-type) for which H3=0, and transverse 
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electric (TE or H-type) for which E;=0t. For the present 
case, these TM and TE wave functions together constitute a 
complete orthonormal solution set from which arbitrary EM 


disturbances in a waveguide can be fully described. 


Hence, a substantial simplification of the solution 
procedure can be achieved by converting the vector 
Eqs.(2.14,2.15) to equivalent scalar equations. TITS 
conversion is accomplished by separating the x;, xX2 
dependent parts of the wave functions (2.18,2.19) into axial 
(E; and H;3) and transverse (E,,2 and H,,2) components. It 


LSQathen pSubticient. to solve 


Vikisiet ysis es 10 (2 On) 

VG ay Hae 0 C22 2u7 
where y?=kg-k3, subject to the boundary conditions 

E3;|Si,2 = 0 (225223) 

0H; = 0 (2233) 

on 54 2 


respectively, since the transverse field components can be 
obtained directly from the axial components according to the 


relations 


tA third species—the transverse electromagnetic (TEM) 
wave—should be considered in a general discussion of guided 
waves. However, the TEM mode cannot exist inside a single, 
hollow duct of infinite conductivity (i.e. the present case) 
and is therefore neglected here. 
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H, = +1Ko 25' (a3 x V,E3) (2Qe24)) 
7% 

E, = + k3 oma a salx H,) (2225) 
Ko 

for TM waves, and 

E, = Fiko Zo (a3 xX VeH3) (2726) 
ape 

H, = a ks Goutal » 4 Es) C221) 
Ko 


for TE waves, where Zo = Yuo/eo, and the plus/minus’ signs 


again indicate forward/backward propagating waves. 


The Eqs.(2.20,2.21), together with the boundary 
conditions (2.22,2.23) respectively, specify eigenvalue 
problems in two dimensions. In general, there will be two 
doubly infinite sets of eigenvalues—one for the TM waves 
and one for the TE waves—that emerge, and the corresponding 
eigenfunctions describe distinct field patterns or 
eigenmodes (or simply "modes") which propagate in the guide. 
However, due LO the mathematical Symmetry of the 
eigenfunctions with respect to differentiation, the 


resulting equation for the eigenvalues is 


Vom=tké sekgee spin? Satp oan: , ('2722'8/) 
L? L3 


where p;,P2 are nonnegative integers, regardless of whether 


T or TE modes are_ sought. Although, the particular 
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combinations of p,,pz which yield nontrivial solutions may 
be different for TM and TE modes. Brome BOs (25.20, Camels 
easily seen that the wave number k; is real only if the 


condition 


Hoe owWrm > Dims + ps ms (2562 9) 
Lj L3 


iSsecatiSived maVvalvueseot p,,p2. thatedomnotesatistyakoq.(2429) 
describe nonpropagating or evanescent modes. Eq.(2.29) thus 
defines the low-frequency limit or cutoff frequency of 
propagating modes. Also note that there is no upper limit 
imposed on propagating waveguide modes. Cavity modes do not 


share this characteristic. 


2.2 Rectangular Cavity Modes 


The formulation given above for waveguide modes can be 
extended to cavity modes by considering plane short circuits 
placed transverse to the x3 axis at x3;=0 and x3=L3. Then 
the domain of interest iS completely bounded by perfect 
conductor and the assumed harmonic x3 dependence is no 
longer valid. The correct form of the x3; dependence is in 
fact the same as that for a standing wave, as would be 
expected for the fields in a waveguide between perfectly 
reflecting planes. Also, an additional set of boundary 
conditions must now be applied at x3=0,L;3. These are 
essentially the same as Eqs.(2.22,2.23), but with E; and H; 


interchanged, ViZ. 
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OE; = 0 (293.0) 
on $3 
H3;|s3 = 0 (243 1) 


which are satisfied for all x; only if the wave number is 


given by 
3 
which is now discrete rather than quasi-continuous as in the 


waveguide case. Also, the relationships between the axial 


and transverse field components are now given [51] by 


Hemaetl Kondo (COS™ Ds UX suilasme VTE) (2h33) 
2 is 
E,. = +p37 sin 3TX3 V.B3 (2.34) 
tia fe L3 


for TM waves, and 


E, = Fiko Zo sin 3TX3 (a3 x V.Hs) Cy SS) 
ve L3 

H, = +37 cos 37X% 3 V.Hs (2.36) 
Lye L3 


for TE waves. 


The cavity case is thus an eigenvalue problem in three 
dimensions and the discrete field patterns (or modes) are 
Stationary and characterized by triply infinite sets of line 
spectraw) Substituting Bq. (2.32) into Eq.(2.28) yields 


ii? jebits, 1 epi ay jake 2rous) 
L } L3 L3 


tates: od pierre ceth conn) 
qe[rez wobee ug", 4278 Hinsonmpa. bialt 
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or, Since kg=€ouow?, 


Uo€oWh = pin? + pin? + pin? C253 G.) 
Lj L3 L3 


which again applies to both TM and TE modes. Note that 
Eq.(2.38) is an equality (compared with Eq.(2.29) which is 
an inequality). The subscript p has been added to the 
angular frequency in order to emphasize that the spectra are 
now discrete. The@ resonant frequencies dW eeare hence 


appropriately called eigenfrequencies. 


The foregoing does not change substantially if the 
Cavity iS completely filled with a perfect (i.e. linear, 
nondispersive) dielectric material. Essentially, the 
QuantitiesS €>» and uo would have to be replaced by the actual 
e and wu which characterize the particular dielectric under 
consideration. However, if the dielectric is not perfect, 
or does not completely fill the cavity, then the boundary 
value problem is changed significantly and the analysis 


presented here is no longer adequate. 


2.3 Rectangular Cavity Field Components 


The eigenfunctions are given (for several geometries) 
in most texts but in order to facilitate later discussions, 
it is useful to list the rectangular cavity eigenfunctions 
in a form consistent with the notation used here. Hence, 


following the procedure outlined in the two previous 
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sections, the actual field components (eigenfunctions) for a 


rectangular cavity are: 


E; = walt: Pim P37 COS 17X 4 sin D27X2 sin D37X3 (23 9D) 
Yiu Dae L's L, Le L; 

E2 = -C p2m p37 sin pimx, COS p2mx2 Sin p37x; (2.40) 
ove L2 L3 L, L2 L3 

E3 = Cy san DPi7X, sin P27™X2 COS 3X3 G2raan) 

L, L2 L3 

H; = 1WE oC p27 sin 17™X1 COS P27TX2 COS P37X3 G24) 
Yicuelts L, L2 L3 

H2 = -iweoC pit COS pix, SiN P2TX2 COS P37X3 C243) 
ope he L, L2 L3 


for TM waves, and 


H, = -D pi® pst Sin pimx, COS p2mxX2 COS P3mXx; (2744) 
Ya L, L3 L, L2 L3 

H2 = -D p2m P3m% COS Pix, SiN P2TX2 COS P3773 (2225)) 
yo bee eLs L, L2 L3 

H; = D coS p,mX; COS P2TX2 SIN P3mX;3 (2554 6:) 

Ly L2 L; 

BE; = 1WUoD p27 Eos 1UX4 sin 27X2 sin P37X3 (20.407) 
gt er aby L, L2 L3 

E> = -1wuoD D,7 sin 17X31 COS TX 2 sin 3TX3 C25) 
gfe nie L, L2 L3 


for TE waves; where C and D are arbitrary constants, and 


harmonic time dependence is understood. 


2.4 Eigenmode Distributions 


Any finite number of the eigenfrequencies defined by 


Eq.(2.38) can be calculated and arranged in ascending order, 
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where w is some arbitrary upper-limiting frequency to be 
considered. The total number of eigenfrequencies not 
exceeding w is then 


N = Ds 1 eqreiey) 


Wp Sw 


where each term in the series represents one eigenfrequency 
in the sequence (2.49). And the toca number of 
eligenfrequencies within the bandwidth 4w is 


Déw = Div ise 2 1 Waoshie) 


Wp swt dw Wp Sw 


where 6 is a (finite) real constant. Note that 
Egs.(2.50,2.51) can be viewed as implicit "functions" of 
frequency Since. 1S) arbitrary sand scan therefore begua 
continuous variable. In this sense, Eq.(2.50) Serves as a 
mathematical definitiont of a spectral distribution of the 
eigenfrequencies. And Eq.(2.51) serves as a definition for 
the spectral density, of . eigenfrequencies, . per Sw of 


PAST Only elite rent derail t1On = 15.) USC somes Ore Lie 
physics literature, viz. 


«ReaD. seine, Dowblie7s 


Wp sw Wp =W 


but this work is mostly concerned with the asymptotic value 
of N in the high frequency limit, in which case the above 
definition and Eq.(2.50) are practically equivalent. 
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bandwidth. Notice that Eq.(2.51) closely resembles the 
definition of the derivative of the "function" N. In fact, 
if 6wx0, then D represents a sum of Dirac delta functions 
5(w-wp) which would be the result obtained from the 


differentiation of the step function described by Eq.(2.50). 


Since Eq.(2.38) holds equally well for TM and TE modes, 
it would be reasonable to define NoweNte and Dry,Dre using 
Eqs.(2.50,2.51) respectively, where only the actual w,) are 
different. Hence, the total number of electromagnetic modes 


not exceeding w would be 
Nem = Nqm* Noe Nn) 
and the number of modes within the bandwidth 6w would be 
Deyow = Drow + Dy dw C2553) 


where Ny,, and N-; are the partial distributions of TM and TE 
modes respectively, and Dy,, and Dyp are the partial spectral 


densities of TM and TE modes respectively. 


It is extremely important to realize that Ny, and Nye 
are not, in general, equal over any given range of 
frequencies, and likewise for Dy, and Dye. This is due to 
the fact that some combinations of pi,p2,p3 uniquely 
identify certain modes as either TM or TE. That is, the 
BolUurtOnmespace Of phd s (2050) ei Seurpalli toned eintoeuhres 


unequal sets of triplets (pip2p3;) which correspond to TM, TE 
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Figure 2.1 Venn diagram representation of the _ solution 
Space of Eq.(2.38) for finite bandwidth. The numbers in 
parentheses indicate the multiplicities with which each 
member occurs in the distributions described by 


Pass(2.52,2.53). 
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and longitudinal (or forbidden) modes. Although the TM and 
TE sets do overlap, there are members which are unique to 
each set. This partitioning of the solution space of 


Eq. (2.38) is depicted schematically in Figure 2.1. 


In order to calculate the mode distributions correctly, 
the members of each of the sets indicated in Figure 2.1 must 
all be known (or at least knowable), which in turn’ requires 
consideration of the boundary conditions (2.22,2.23) and 
(2-30) h2053 1) These mathematical expressions can be 
translated into restrictions on the quantum numbers that 
COLLeSDONGsNLO =NOntLri Vial —SOlULTONS= at Om mEOS.\ 26 205s cs 
These restrictions will be described in detail in the 
following chapter where they are formalized into a set of 
Criteria for properly discriminating between TM and TE 


modes. 
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CHAPTER 3. NUMERICAL EIGENMODE DISTRIBUTIONS 


Several previous attempts to enumerate the eigenmodes 
in empty rectangular cavities by numerical computation have 
produced incorrect results. The purpose of this chapter is 
to correct these defective algorithms. A brief review of 
the physical principles involved in electromagnetic mode 
counting is given in the first section. A particular mode 
counting algorithm which yields the correct number of modes 
in a rectangular cavity is presented in the second section. 
Three different implementations of the algorithm are briefly 
described in the third section, followed in the next two 
sections by some important numerical results obtained from 
two of these programs. A discussion of the significance of 


these results is given in the last section. 


3.1 Numerical Mode Counting 


The spectral distribution of modes for a particular 
cavity may be presented graphically as a plot of the total 
number of modes with resonant frequencies not exceeding some 
upper frequency w as a "function" of w. Similarly, the 
spectral density of modes can be presented as a plot of the 
number of modes with resonant frequencies between some lower 
frequency w and wtéw as a "function" of w. In either case, 
the computational problem essentially reduces to counting 


the number of allowable modes (eigenmodes) within a 
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specified frequency range. For rectangular cavities, these 
mode countS can easily be obtained using exhaustive 
substitution algorithms of the type described in the 
following sections. In=™= fact; @the®rresults# of Similar 
computations have been reported in the form of tables 
[42,43,48] intended to aid designers in the choice of 
optimum cavity dimensions. However, the algorithms used by 
these authors do not account for the physical differences 
between TM and TE modes, and therefore do not yield the 


correct total number of modes. 


The intrinsic eigenfrequencies of an ideal rectangular 
Cavity are defined by Eq.(2.38) which, in the engineering 


literature, is usually written 


2 th Getter Xo jepke’ “2 yseie- G3 Fat) 
4L? AE,2 4L3 


where £ is the more commonly used Hertzian frequency, c is 
the velocity of light, and the quantum numbers p;,,p2,p; 
together comprise the so-called mode number (pipzp;). 1 
there were no restrictions on the mode number (except that 
it be nonnegative), then Eq.(3.1) could be used to generate 
the aforementioned plots by simply counting all of the 
possible mode numbers which define resonant frequencies 
within the appropriate’ range. Thi.SaielS MeL ieecac CrammEEne 


approach used to compile the tables in [42,43]. 
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As explained in the previous chapter, Eq.(3.1) results 
from the solution of the Helmholtz equation regardless of 
whether TM or TE modes are_ sought. Hence, in the 
hypothetical treatment described above, one might assumet 
that the total number of modes is just twice that determined 
for either species alone. Such a treatment presupposes that 
TM and TE modes have identical partial distributions, which 
theysdo not... That »is, the TM wave. function® Satisfies ) the 
Dirichlet ‘condition! on' those portions s,of ther sboundary 
surface where the TE wave function satisfies the Neumann 
condition, and vice versa. As a result, if any one of the 
quantum numbers is zero, then one or the other wave type can 
Exnistwambutmenct aboth< And if any two or three quantum 


numbers are zero, then neither wave type can exist. 


Notwithstanding the above, certain mode numbers are 
ambiguous in that they can represent both a TM and a TE 
mode. Thus, any algorithm which is intended to count 
electromagnetic modes must include a facility for assigning 
multiplicities inmiorderwhtom accountgefors thismetypemsot 
degeneracy. That is, 1f none of the quantum numbers are 
zero, then the mode number represents both a T and a TE 
mode, and therefore can be said to have a multiplicity of 


two. Whereas, a Single zero occuring in any position in the 


+The instructions given in [42,43] are not explicit, but it 
is obvious from the context and the actual results that only 
one type has been considered, and no distinction between 
types is made clear in the discussion. 
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mode number will uniquely identify the mode as either TM or 
TE, which therefore has a multiplicity of one. The 
remaining case of two (or three) zeros in the mode number 
leads to the trivial solution where all of the 
eigenfunctions vanish, and so these modes are said to have a 
multiplicity of zero. These multiplicities are thus the 


same as those indicated in the diagram of Figure 2.1. 


Due to the simplicity of the eigenvalue problem in 
rectangular geometry, the numerical analysis is 
Straightforward (perhaps even elementary) once the physics 
is understood. Still, in order to establish a basis for 
comparing theory and experiment, it is worthwhile to examine 
a particular algorithm based on the above physical 


conditions. 


3.2 A Mode Counting Algorithm 


The algorithm consists of a simple ide ange 
procedure for determining which of the possible mode numbers 
yield resonant frequencies in the desired range. And a 
"Sorting" procedure for determining which of the filtered 
mode numbers correspond to valid cavity modes, and whether 
they can be TM, TE or both. The filtering procedure is 
essentially just an application of Eq.(3.1) and does not 
merit further consideration here. Hence, only the sorting 


procedure will be described in detail. 
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A thorough sorting procedure (i.e. one which can 
discriminate between TM and TE modes) requires a slight 
modification to the criteria briefly described in the last 
section. The actual criteria used in the present algorithm 
for sorting modes according to mode number are the 


following: 


(a) If more than one of the p,,p2,p3 are zero, then the 
mode number represents neither a TM nor a TE mode (i.e. a 


longitudinal, or forbidden, mode). 


(b) If only p3; is zero, then the mode number represents 


a TM mode only. 


(c)Meert ) either py, or p, (but@not p, iis) zero.) then the 


mode number represents a TE mode only. 


(d) eb fsanonegeofMathGapgsp2yps aregzero, then the mode 


number represents both a TM and a TE mode. 


The modification indicated by (b) and (c) follows 
Girectly from the complementarity of the boundary 


conditions. 


Let P,,P:,P3 be Boolean (input) variables which are 
true if p;,P2,P3 respectively are nonzero, and false if they 


are zero. Also let E, H, B be Boolean (output) variables 
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which are true if a given set of input variables describes a 
T (or E-type) mode, a TE (or H-type mode) or both 
respectively, and false otherwise. These definitions 
together with the criteria (a)-(d) can then be used to 
construct Table 3.1 which comprises all of the possible 
combinations of P,,P2,P3 and the corresponding binary values 


ope fe 456 Ae 


In order to obtain an efficient algorithm, the minterm 
expansions for E, H and B are derived from the Karnaugh maps 


Sa easmsSnowneine Tablewsa2s) V 172% 


E = P,P2 Sha Zh) 
H = P,P; + P2P,\= (P; + PoP; C3553) 
B = P,P2P;3 (3.4) 


which can then be used to test and validate (or invalidate) 
the mode number (piP2P3). Similarly, a test for 


nonlongitudinal modes can be derived, viz. 
N = P,P; + P2P; + P,P2 (gk 


where the definition of N is analogous to those given for E, 


Hees eeNOtesthateNenls sJ)UStmnthHbas would spe wexpected, . 


+Perhaps all of the Eqs.(3.2-3.5) could have been written by 
inspection but, inasmuch as Karnaugh maps are a _ proven 
mathematical construct, the procedure described can, given 
Table 3.1, serve as a proof of these equations. 


7 2 


At 


. 
precipi hides abs atin ake “— a 
eiMiveug! nets! dks ese lem Savas tact 
soviev ycanid aetbaaqessi6s “ou CTH FT 26 2m | 


~ 
4 


a : 
wh 
sheet ie G ; 


* 


_ 
Mrernie ods aitsixgote ie 8 a rietds ct t6bae. a 
. 

Pi 


ss 


gent eplawis? af? tort havete alia roe a 
yy 7é “lav st avots am t i} 


ah. =) sti4 a g 
t€,@) 3.9 2493) = 63.3 2 = 8 z 
$5? 63e9-S = @ 


(eystiievuk.s0) otabifsy Pode. Jeaes oF beeu $4-pnytt? 2a9 e} iw 
743 Jaa3 & ¢fzstlmi2 »(eoe0,a) \eceun oboe 4 {3 


2hy> ba¥iasb" a4 fixe! asen fanthes ignetieae 
—s : 7 : : 

i : : . 

(9.4) 8% (957 * pet ow 


: 9 as 8 =. , ; sé = aan 
8 10] asvip seadd oP -ayopufans ai “do. udternidee ate eat 


tosioaqey ad Bluow ge Hee zed) 2i i dda Sa G 


31 


Table 3.1 Truth table for the Boolean output variables E,H,B,N 
which are true if the corresponding input combination describes a 
TM mode, a TE mode, or both, or a non-longitudinal mode 
respectively. The input variables Pi rPorP are true if the 


quantum numbers Pi +Po/P3 respectively are non-zero. 


“This column describes the algorithm used (implicitly) in [48] 
which yields different, though still incorrect, results from those 
reported in [42,43]. 


Din Slight simplification may seem possible if these are designated 
as "don't care" states, but a more complicated filtering procedure 
would then be required in order to eliminate these cases. 
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Table 3.2 Karnaugh maps for the Boolean output variables 
E,H,B,N. The resulting logic equations are written in minimum 


sum-of-product form as described in [54]. 
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Eqs.(3.4,3.5) can provide additional information if 
required, but Eqs.(3.2,3.3) are the essential core of the 
sorting procedure. The actual number of TM and TE modes, 
and therefore the total number of EM modes, can thus’ be 
counted during or after sorting. Ai low Vcharteewhich 
describes a likely implementation (and one used by the 
author) is presented in Figure 3.1. This algorithm was used 
to generate numerical data for a particular laboratory 
cavity for which experimental data could be obtained for 
comparison. These results are compared graphically in the 
next Chapter, for now let it suffice to say that the results 


are in close agreement. 


3.3 Computer Aided Cavity Design 


As expected, the above algorithm is not very 
complicated, and more efficient algorithms no doubt exist, 
at least for some special casest. But for practical cavity 
sizes and bandwidths, used in microwave heating 
applications, the computation time for any one cavity is 
reasonable (i.e. less than one minute, using interpreted 
BASIC on an 8 bit microcomputer; and less than one_ second, 
using Pascal or FORTRAN on a 32 bit mainframe). This 
algorithm has in fact been implemented, with satisfactory 
results, in a trilogy of interactive design/analysis 
+For cubic cavities for example, which are inherently the 


most mode-dense, a much more efficient algorithm is 
described in [27]. 


toasts mugen A 
sits td beau SAS Gite) 
fers acy ada tnorte. abit 
qrodetedat seluae ree 
102- bed lsdee os athe: éynb. * $eerant2s95* 


chbsetaaars fowi 


tordw 
of> «i ech sebicdge oe aadmon S30 erivees eteeT 


2jizea: add zens. Hee, Sf estPide 4i Jel Wom*302 


itengeipe ecois ab oe 


npiasG ys:vad WektA tetped Bb 


y?sv me 2 isbn ae aveve ons 


fates shaae on 2ad P=tile augico: tts sicm af 


7 he 4023 TR SaGS5 


Eau 


yrives. fap tiseck 
dees Ao es 
at acevan ano, jou 
yal -piges ra ant uted a5 
ae ng 


8h *bivboed 


qPeRl? nos te4URReD aa? 


“AG THON 30 


tae tes so or 
yisail = @ 


ie giag!t af feaqeetg st una ‘5. 


Pe 
he Pe 2 arab  feditsnun staeanep c P 


\Sercuges GA 


($ipage2 ence 207 Jasel 
S46. 


eek} 

Bhs (iszuytosorsia 2id 6 as fo 21 
: _— 

esd shea ft ena nase 


= 


ele 


13 estes 


-cene 7 


,Tsiget pe - 


7 > f 


nna 


34 


FILTERING PROCEDURE 


EQUATION (3.1) 


———— _ — ee eee eee =e es 
Seeeeee= — wee eee | He ee me me em ee ee eee ee ee we ee ew ee me Oem ee eee 


EQUATION (3.2) INCREMENT T™ 


MODE COUNTER 


EQUATION (3.3) INCREMENT TE 


MODE COUNTER 


THE TOTAL NUMBER OF 
MODES IS THE SUM OF 
THE T AND TE COUNTERS 


Figure 3.1 Flow chart representation of a mode counting 
algorithm for rectangular cavities. The applicable 
equations from the text are indicated beside the 


decision-making blocks. 
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programs used by the author. 


In the program NEWCAV, the algorithm is used to find 
the most mode-dense cavities (over a given range of 
frequencies) in a set of cavities defined by varying, ina 
prescribed manner, a set of target dimensions. The program 
prints a listing of the eigenfrequencies, along with plots 
of the TM and TE line spectra,. for each of the cavities 


found. 


In the program MODIST, the algorithm is used (along 
with certain of the analytical formulae presented in Chapter 
5) to generate mode distribution data with any required 
degree of resolution. The output data is written, in. 
tabular form, to a file which can be viewed as is, or used 


as input to a plotting program. 


In the program CARETAB, the algorithm is used to 
produce mode count reference tables for a set of cavities 
defined in a similar manner to that used in NEWCAV. The 


output of this program is described in Section 3.5. 


Complete listings of each of the above programs and 


descriptions of the input/output are given in Appendix A. 


A distinct limitation of these and similar CAD programs 


used in engineering is that the presence of the coupling 
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structure(s) and mode stirrer perturb the cavity enough that 
EQuiseal) mUSenomerongermmuselful A considerable advantage 
would be gained if a numerical model could be devised which 
is sufficiently general to handle multiple antennae and 
arbitrary metallic structures within the interior of the 
cavity. It is beyond the scope of this thesis to consider 
this —topic in “detail; but some further thoughts on a 


possible approach are given in Appendix C. 


3.4 Numerical Case Studies 


The program MODIST described above can be used in 
conjunction with a graphics program to generate actual plots 
of the exact spectral distribution or spectral density of 
modes for any hypothetical cavity, for any finite range of 
frequencies. It is interesting to examine the results of 
such numerical case studies and some particularly 


instructive examples are presented here. 


The ‘cavity depicted in Figures 3.2,4,678) 1S cubic 
(L,=L2=L3=25 cm) and therefore completely symmetric under an 
interchange of coordinate axes. This case gives rise to the 
maximum number of degenerate modes, and hence to _ the 
greatest step heights in the exact spectral distribution 
(see Figure 3.2). Thisveine tUrn@edeads @toueethetalargest 
fluctuations in the spectral density plots (see Figures 
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Fs U iT Cameedie The theoretical spectral distribution of 
electromagnetic modes in an empty cubic cavity. The 
dimensions of the cavity are L,=L2=L3=25 cm which give rise 


to a cutoff frequency w,/27=848.5 MHz. 
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Figure sae The theoretical spectral distribution of 
electromagnetic modes in an empty rectangular cavity. The 
dimensions of the cavity are L,=25 cm, L2=35 cm, L3;=45 cm 


which give rise to a cutoff frequency w,/2m7=542.9 MHz. 
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Figure 3.4 The theoretical spectral density of 
electromagnetic modes in an empty cubic cavity for a 
bandwidth of 100 MHz. The dimensions of the cavity are 
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Figure eh) The theoretical spectral density of 
electromagnetic modes in an empty rectangular cavity for a 
bandwidthesoLe 100 MMHz. The dimensions of the cavity are 


Dye25 em, G2=35ecmn,eb5,=45 cm as in Figure 3.3, 
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Figure 37.0 The theoretical spectral density of 
electromagnetic “modes in an empty (nearly) cubic cavity for 
eebandwidtheor  1005MEZz.= “The dimensions ot “the =cavity: ‘are 


L,=26 cm, L2,=L3;=25 cm. 
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Figure Be The theoretical spectral density of 
electromagnetic modes in an empty rectangular cavity for a 


Dandwadthy got 901008 MHZ: The dimensions of the cavity are 


L,=26 cm, L2=35 cm, L3=45 cm. 
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Figure 358 The theoretical spectral density of 
electromagnetic modes in an empty cubic cavity for a 
bandwidth off 50 MHz. Thera dimensionss,of sthesscavitv, are 


L,=L2,=L3;=25 cm as in Figure 3.2. 
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Figure sae) The theoretical spectral density of 
electromagnetic modes in an empty rectangular cavity for a 
Dandwidthwe of.) 50, .MHz,. The dimensions of the cavity are 
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Symmetric (L,=25 cm, L2=35 cm, L3;=45 cm) but the step 
heights in the spectral distribution (Figure 3.3) and the 
fluctuations in the spectral density (Figures 3.5,7,9) are 
substantially reduced. These figures clearly show the 
marked dependence of the mode distributions on the degree of 
symmetry of the Cavity under consideration. Notice 
particularly the pronounced difference between Figures 3.4 
and 3.6, and between Figures 3.5 and 3.7, where the only 
change in the parameters is that the L, dimension has_ been 


perturbed by 1 cm. 


Figures 3.8 and 3.9 show another strong dependence 
(that is, on variations in bandwidth) which affects only the 
spectral densities. A visual comparison of Figures 3.8 and 
3.9 with Figures 3.4 and 3.5 respectively, indicates the 
sensitivity of the fluctuations in the spectral density toa 


change in bandwidth (in this case, from 100 to 50 MHz). 


3.5 Revised Cavity Design Tables 


For the sake of comparison, example mode count data,’ 
obtained from the program CARETAB, is given in Table B.1 of 
Appendix B, in a form similar to the tables given in 
[42,43]. All entries pertain to the number of eigenvalues 
of the specified type within the (Hertzian) frequency range 
2425 s —£ < 2475 MHz. The column titled TS is the total 


number of mathematical solutions to Eq.(3.1), including 
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double-zero (longitudinal) modes. These are the values 
given in [42,43]. The column titled NL is the total number 
of nonlongitudinal modes (i.e. the same as TS but with the 
double-zero modes extracted). This is the method used in 
[48]. The remaining columns titled TE, TM and EM are the 
partial and total mode counts as determined by the present 


algorithm. 


A Gsinilar §seterofigdatay! foraithe frequency range 
S20Cetst Vi veser 930MeMHZ areugivensin&éTablegBt2 offAppendix Bs 
Note that, due to the nature of the exact spectral 
distribution of modes, it is not possible to interpolate 


between cavity sizes given in Tables B.1 and B.2. 


3.6 Discussion 


It is not clear whether or not the authors of [42,43] 
intended users to double the values given in the tables. If 
they are doubled, then these results consistently 
overestimate the total number of modes. But if they are not 
doubled, then these tables would consistently underestimate 
the total number of modes, which may explain why the 
coupling to small cavities is often better than expected. 
In either case, the discrepancies are quite large. Note 
however, that these discrepancies may not have been possible 
to detect experimentally at the time of [42,43]. Since 


solid state microwave applicators will almost certainly 


7 


ae 3 
. ‘asin ine tives eds vt It.b5 ales mufos st TOD, eT aE 


seats net Som (fen ttt vet a1 
i ot 


baal Ps 


7s - 
a 


sad d3tw tut oy ¢4 she its er) gabon. wire: jinen 20 


yi - 7 i J du 
al bass Bofgem off 84 eet >the e15H7I49 © om = _ 


;' bet ois sningd abe on tnken or. ry 
: _ 


= 


ant? e716. 42 ore MY .5} 


pia ti Wed acen Tatod tae 
oe .— =e an 4 ramerints 4 se 3 hehe ale pe A305 {ajod r 
#Pses- ‘] af 1.» a to 1) ie - { 
a vines ize aq? “107%, (&350 isk “TSi se 


~~ - o> ~~) E ¢ 
ta.200¢2 25625 Sri? 70 ‘Tu or 3.1 wu wo 4 
gvatoteasnt of Sigideta-donves; 7:  eshom x 
- ryt b : 2 
<€.8/ Bag %o"-adgileaT ar nsvie sill vitived neat 
= ~~ “ 
7 | j.BeuGg 
yy : ~T <s-o7 PP 
fea Shi = Yo esedsus of. 10n 3c t8hIah¥ T9sA9 ef 
2t. ,sstded sf2 of cavin wsulav sic efdves cy asaae & 
SS oes ay rt , 7 : 
% * or. ta = 
(1431S TEI EAOT 0, gta et SSab,? igs _ \ coe sete 
: @ ‘ - a : 7 — 
- - Py ~ 
foe) S76 y9ns <=! ue ,2sooq 356 Tasnen -Lador ant stomi jas ovo 


sinwiztveony yisaeda aie Btu Saaiteiatt -esans vad) bs dgvob 
- : ez - as * : ava 
edn did. Aieteue [yen fpiev ,eahon tc — lasod off, 
Rint oo 2i asi zivas: fta mes 12 £que: 
7 


wep is 
re _ 


47 


operate at lower frequencies initially (e.g. the 915 MHz ISM 
band), and with multiple sources, the coupling problems of 


small cavities are becoming increasingly more important. 


It might seem that the existence of quasi-random 
fluctuations in the spectral distribution demonstrates the 
necessity for designers to have available extensive mode 
distribution data (or™at™ least "a *facility™ to obtain such 
data) in order to optimize cavity dimensions. However, the 
sensitivity of these fluctuations to slight variations in 
shape and bandwidth raises some uncertainties regarding the 
interpretation of data obtained for ideal empty cavities. 
For example, construction tolerances and dynamic magnetron 
bandwidths cause unpredictable changes in the spectral 


density of modes in "real" cavities. 


Furthermore, the resonances are quite sharp, even in a 
(lightly) loaded cavity, and sometimes widely spaced. 
Hence, in order to achieve satisfactory uniformity and 
efficiency in microwave applicators, mode stirrers and/or 
turn tables are usually employed. These drastically perturb 
the shape and symmetry of the cavity, which substantially 
changes the actual mode distribution. Thus, even if 
intrinsic mode distributions are accurately obtainable from 
numerical computations, one cannot necessarily expect to 
improve on empirical cavity designs. Some of these aspects 


are examined experimentally in the next chapter. 
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CHAPTER 4. EXPERIMENTAL EIGENMODE DISTRIBUTIONS 


The algorithm described in the preceeding chapter is 
based on principles which are strictly valid only for ideal 
empty cavities. The purpose of this chapter is to present 
the results of experimental measurements of the spectral 
distribution of modes for a "real" cavity at microwave 
frequencies. The laboratory cavity, the measurement system 
and the procedure are described in the first section. 
Experimentally obtained mode distributions for the empty 
Cavity are compared in the next section with those obtained 
by numerical computation using the algorithm of Chapter 3. 
Additional measurements on the same cavity, equipped with a 
mode stirrer, are presented in the third section. These 


results are;discussed,inithe last section. 


4.1 Materials and Methods 


The laboratory cavity is constructed of aluminum with 
InteruoOnendimensl1 ONS eel s=com CN aeus=00 CMy4lb5=52ecmenandecan 
accommodate tunable probe type coupling structures on three 
orthogonal’ sides. Access to the interior of the cavity is 
provided via a securely fitting side door (23 cm x 52 cm), 
or through several small "ports" (1 cm dia.) which can be 
used to introduce perturbing instruments or additional 
probes. An aluminum mode stirrer, of the type used in 


Litton Menumastero microwave ovens, could be installed in 
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the top wall and rotated manually from outside of the 
cavity. An aluminum bushing was used in order to maintain 
electrical contact between the mode stirrer and the cavity 


wall. 


The distributions were measured by counting the number 
of minima in return loss vs. frequency traces obtained using 
a coaxial frequency-domain reflectometer system. A 
swept-frequency source (HP 8620C, HP 8621A, HP 8632A) 
capable of a continuous sweep from 0.1 to 4.0 GHz was used 
to excite the cavity. Incident and reflected power signals 
were sampled using -20 dB directional couplers, and the 
return loss in dB was displayed on a swept amplitude 
analyzer (HP 8755A). No attempt was made to match the 
impedance of the input circuit to the cavity in order to 
keep the loading effect of the external circuit as small as 
possible. A schematic diagram of the complete measurement 


system is shown in Figure 4.1. 


Degenerate modes were resolved using perturbation 
methods [55]. In many cases the perturbing effect of the 
finiresconductivity Of tne walls and™ the “presence mot sethe 
coupling probet were sufficient to split degenerate modes 


just enough to make identification possible. From the case 


+At least two probes (in at least two positions) were used 
in any given measurement run in order to excite all possible 
modes, but only one probe was excited at a time and the 
unused probes were removed from the cavity. 
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SWEPT AMPLITUDE ANALYZER SWEEP OSCILLATOR / RF SECTION 


z-axis in 


=. 


_—————————— 


external in Sweep out 


modulator drive 


DUAL DIRECTIONAL 
COUPLER 


MODE STIRER 


Figure 4.1 Schematic diagram of the measurement system used 
to obtain experimental mode distribution data. The inset 


shows the location of the mode stirer used for the 


measurements reported in section 4.3. 
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studies discussed in the previous chapter, one would expect 
a higher sensitivity to probe perturbations in a more 
symmetric cavity. Some abbreviated measurements (not 
graphed) on a slightly more symmetric cavity 
(28 cm x 52 cm x 52 cm) did support this, as extremely short 
probes (<5 cm) were required in order to observe degenerate 
modes. However, the coupling from these short probes was 
extremely weak which made some degenerate modes very 
difficult to discern. Because the probes could be longer 
(<10 cm), this problem was less serious for the cavity used 
in the measurements reported here; although distinguishing 
weakly-coupled degenerate modes was still the probable cause 


of most of the error due to "missed modes". 


Other sources of "miscount" error include the possible 
Gountinge lof gartifactSyeaing thes returnjseioss trace and 
recounting of modes that are excited by different probes. 
The artifacts could be identified in most cases by 
repetitive measurements and perturbation methods. The 
chances of recounting some modes could be lessened by 
exciting all probes simultaneously. However, cross-coupling 
between probes then introduces another source of error. 
Some (qualitative) observations of these effects are 


reported in Appendix C. 


Errors in determining the resonant frequencies of 


nondegenerate modes and degenerate modes that required 
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manual intervention to split, were mainly due to the limited 
resolution of the scale (on the sweeper main-frame) used to 
obtain the frequency Gata. However, the resonant 
frequencies of modes that were split due to the perturbing 
effect of cavity walls and coupling probes could only be 


estimated. Fortunately, this type of error is small. 


4.2 Results for the Empty Cavity 


The experimentally measured mode distributions for the 
empty cavity are plotted in Figures 4.2 and 4.3 as _ the 
curves N™ and D™. The abscissae are all normalized to the 
cutoff frequency (i.e. the lowest eigenfrequency) w, of the 
Cavity. Superimposed on Figures 4.2 and 4.3 are the 
(numerically determined) theoretical distributions N‘t and D* 
respectively obtained using the algorithm presented in 


previous chapter. 


The frequency resolution of the computed distributions 
was chosen to correspond with the resolution of the 
experimental measurements, which was estimated to be 10 MHz 
or approximately 0.023 on the normalized scale used in 
Figuresmé2eanda4,.d.——eWibhinethiseresolutionjetheresisecilose 
agreement between the computed and corresponding measured 
data. The deviation at higher frequencies is attributed to 
the experimental errors mentioned above. Note that the 


visual counting of modes becomes increasingly more difficult 
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Figure 4.2 Comparison of the theoretical (Nt) and measured 
(N") spectral distribution of modes in an empty rectangular 
Cavity. The interior dimensions of the cavity are L,=23 cm, 


De=ouecm, Gs-o2 cm which Qive rise tow a "cCuvot rt *trequency 


wo /2m=416.2 MHz. 
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Figure 4,3 Comparison of the theoretical (Dt) and measured 
(D™) spectral density of modes in an empty rectangular 
Cavity for a bandwidth of 50 MHz. The interior dimensions 
OtmuneacCavityeare (u7=235g0ct) eeus=o0n—BCM ae. .=o2 scm ema Samah 


Figure 4.2. 
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at higher mode densities and therefore the accuracy of the 
measurements deteriorates rapidly at frequencies beyond the 


range shown in the figures. 
4.3 Results for the Cavity with a Mode Stirrer 


The experimentally obtained spectral distribution of 
modes in the cavity fitted with a mode stirrer is plotted in 
Figure 4.4 as the curve labelled NT. The empty-cavity 


distribution N™ is redrawn on this figure for comparison. 


The shaded portions of the curve indicate the range of 
resonant frequency shift due to changes in the angular 
displacement 6 of the mode stirrer. There is a substantial 
degree of uncertainty in the determination of the endpoints 
of these ranges for some of the higher frequency modes. 
That is, the resonant frequencies of adjacent modes 
sometimes shifted to where they momentarily coincide; 
identification of the emerging modes was then difficult 
since the frequency shifting is not necessarily 
single-valued with respect to the angular position of the 
mode stirrer. Also, several of the modes exhibited drastic 
and abrupt fluctuations in coupling due to changes ok Uy Ge 
fact, many modes which were strongly coupled at some values 


of @ could be made to completely vanish at others. 
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Figure 4.4 Comparison of measured spectral distribution of 
modes in a rectangular cavity with (NJ), and without (N™), a 
mode stirrer. The interior dimensions of the cavity are 


De=csecm ele =50 Cm, Le=o2ecem as ine Figures4. 2. 
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Note that the range of frequency shift was determined 
for each mode (or set of degenerate modes) separately, and 
hence the actual 6 dependence of the ensemble distribution 
can not be determined from Figure 4.4. In this sense, N™ 
does not really depict the spectral distribution curve 
per se, but rather it depicts the envelope of possible 
Spectral distribution curves within the frequency range 
shown. It is therefore pointless to attempt to extract any 
mode density information from Figure 4.4 as interpretation 
Ofpesuchreagn curve” Avouldienot be possible. Similar 
considerations apply to the experimental determination of 


the spectral density. 
4.4 Discussion 


The good agreement between the curves N™ and N* and 
(hence) between the curves D™ and Dt constitutes a 
successful test of the algorithm presented in the previous 
chapter. This result is not Surprising since the algorithm 


was derived from theoretically sound physical principles. 


What is perhaps surprising is that the inclusion of a 
mode stirrer seems to have caused a slight negative shift of 
the entire spectral distribution curve (i.e. downward in 
frequency). Since the mode stirrer is located near the 
walls of the cavity, one might expect that magnetic fields 


would be perturbed more than electric fields which should 
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result in a net positive shift in the resonant frequencies 
ES cult. However, the integral of the stored magnetic field 
energy does not exceed the integral of the stored electric 
field energy over a large enough region of the mode stirrer 
volume to result in a net positive frequency shift. 
Although, a slight positive shift could be obtained for most 


modes, for some angular positions of the mode stirrer. 


The perturbing effect of the mode stirrer also caused a 
more thorough splitting of degenerate modes, thus reducing 
the amount of manual splitting required. Furthermore, the 
mode stirrer appeared to act somewhat like a short-circuited 
"dummy" antenna which received power from the exciting probe 
and reradiated it into modes normally inaccessible to _ the 
exciting probe. Thus, all of the modes could be excited 
from a single probe position which was advantageous from a 
measurement point of view. However, this effect would 
virtually eliminate at least one of the proposed methods 
[48] for reducing the cross-coupling between multiple source 


antennae—that of cross-polarized antennae. 


An effect of the mode stirrer which is not conveyed by 
the curve in Figure 4.4 is the marked reduction in the Q of 
Some of the modes even in the absence of other forms of 
loading. This effect may be attributed, at least partly, to 
the decrease in volume and the increase in surface area of 


the cavity due to the presence of the mode stirrer 
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(i.e. since Q is proportional to the ratio of V/S). Also, 
from an equivalent circuit point of view, the mode stirrer 
contributes a capacitance which is in series with the usual 
model capacitance, and since the Q is proportional to the 


total capacitance, a lower value results. 


The experimental methods described here are not limited 
to measurements on rectangular cavities. In fact, some 
interesting effects have been observed in odd _ shaped 
(nonrectangular) acoustic cavities [55] which may be 
relevant to microwave heating applications of 
electromagnetic cavities. The authors of [55] present 
analyses of the statistical dispersion in the excitation 
amplitudes of neighbouring modes, which suggest that more 
uniform field patterns might be obtained with completely 
asymmetric cavities. Unfortunately, experimental results 
from these cavities are difficult to generalize, as very 
little theoretical analysis is possible. Presently, the 
only analytical aid for determining the mode distributions 
for these (nonrectangular) asymmetric cavities is offered by 
the shape-independent volume term in the asymptotic formulae 


discussed in the following chapter. 
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CHAPTER 5. ANALYTICAL EIGENMODE DISTRIBUTIONS 


Asymptotic expansions of the eigenvaiue distribution 
have been used in physics to approximate the exact spectral 
distribution of modes in several types of resonant cavities. 
Certain versions of these formulae have been incorrectly 
applied to microwave cavity problems in engineering. The 
purpose of this chapter is to correct these errors and to 
comment on the utility of the approximate formulae in 
engineering. A brief introduction describing the errors is 
Givensin = sthe (‘first Section, followed in ‘the next “two 
sections by derivations of the correct electromagnetic 
formulae. In the fourth section, these formulae are 
compared with the numerical and experimental results 
presented in the last two chapters. In the fifth section, 
exact formulae for the spectral density of modes in 
rectangular and cubic cavitieS are given. Inthe last 
Section, all of these formulae are discussed in terms of 
their role in the design and analysis of microwave heating 


cavities. 


5.1 Asymptotic Eigenmode Distributions 


Recall that the time independent behavior of 
oscillating fields in bounded domains is governed by the 


Helmholtz equation, VizZ. 
V7y + k?p = 0 (591) 
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where the wave function w may be either scalar or vector 
depending on the type of wave phenomenon under 
consideration. For example, the scalar form of Eq.(5.1) 
applies to (gas phase) acoustic waves and the vector form 
applies to electromagnetic waves, as discussed in Chapter 2. 


In either caset, the eigenvalues of Eq.(5.1) satisfy 


k? = w? = jhe ee fee ce jaime (52) 
C ae ine Li 


where the valid combinations of p,,p2,p3 depend on the 


boundary conditions. 


The previous two chapters both dealt with methods of 
determining the exact eigenvalue distributions as defined in 
Ghapten., 2. Such methods can only give qualitative 
information about the general behavior of eigenvalue 
spectra. However, analytical functions of frequency can be 
derived which fit the curves of N,D empirically, at least 
for asymptotically large frequencies. These functions have 
been used by physicists, for example, to elucidate certain 
properties of blackbody radiation in finite cavities 
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several forms of asymptotic formulae must exist in order to 


+Note that not all problems involving Eq.(5.1) lead to 
Eq.(5.2). A case in point is that of acoustic modes in 
solid-filled cavities due to a nonisotropic phase velocity. 
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handle different sets of boundary conditions, geometries, 
etc. For example, the leading terms of an asymptotic 
expansion of the eigenvalue distribution valid for the 


Scalamerlorm #ofmtbqe( 5 71)) cin la rectangular parallelepiped 


domain is 
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where V is the volume of the cavity and S is the interior 
surface area. The plus signs in Eq.(5.3) apply to Neumann 
problems and the minus signs apply to Dirichlet problems. 
The constant term may be viewed as being necessary to 
properly Vaccount for the 000 mode, but its contribution to 
the total number of modes is usually minor enough to _ be 
omitted. These appropriate eflorm -OfeeEG.G>-s))ukvallcmmcor 
acoustic modes (in a gas-filled cavity) is 


NERC) ot=, Vil? a, S 
acous eons =C. 167 


ce: (L, #D34+D)) Ww - di (5 34) 
eC 8 


which is the form derived in [40]. Since Eq.(5.2) applies 
to both the scalar and electromagnetic problems, it was 
(apparently) assumed that Eq.(5.4) also applies to the 
electromagnetic problems. The form suggested in [4345 i is 
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where the w dependent terms in Eq.(5.4) have simply been 
doubled to account for both TM and TE mode types—a 


procedure which is correct only for the volume term. The 
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problem being that, while the scalar problems admit 
solutions which are oscillatory in one dimension only 
(i.e. longitudinal), electromagnetic problems do not. 
Furthermore, although the solution spaces of TM and TE modes 


do overlap, there are solutions which are unique to each 


set. 


Although electromagnetic waves are governed by the 
vector form of Eq.(5.1), the (vector) wave functions can be 
obtained in terms of two scalar wave potentials representing 
T (E-type) and TE (H-type) modes’ respectively. This 
formulation enables the TM and TE modes to be dealt with in 
separate (Scalar) problems, each subject to a different, 
albeit complementary, set of mixed boundary conditions. 
This complementarity leads to the cancellation of the 
surface area dependent term in the electromagnetic 
aistribution., That is, the terms proportional “to (os for 
each mode type have equal but opposite signs; and hence, 
when the partial solutions are superimposed, the _ surface 
term vanishes. The first rigorous treatment where TM and TE 
modes are considered separately was given in [27]. The 
important results are derived here by a method similar to 


that used by Maa. 
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5.2 Asymptotic Spectral Distribution of EM Modes 


Eq.(5.2) can be interpreted as describing the boundary 
of an ellipsoid in the space of the quantum numbers 
Pi,P2,P3f. Hence, any point in p-space having integral 
coordinates that lie interior to the ellipsoid will 
represent a mode with eigenfrequency less than w. There 
will be, on the average, one of these characteristic points 
per unit volume in p-Space. The number of TM or TE _ modes 
with resonant frequencies lesseythan. @, atom ame rse 
approximation, is then simply the volume in the first octant 


enclosed by the ellipsoid, viz. 


67 2 3 

where the subscript 0 indicates that Eq.(5.6) is applicable 
to both TM and TE modes. Eq.(5.6) is the familiar Rayleigh 
formula and it (or a multiple of it) is the dominant term 


which describes all eigenvalue spectra in the limit of 


infinitely high frequencies or infinitely large cavities. 


Note that the mode numbers comprised only of nonzero 


quantum numbers have been fully counted in Eq.(5.6). 


tThe hypothetical surface described by Bq etd. 20 is 
ellipsoidal due to the simple rectangular geometry of 
x-space. Clearly, more complicated surfaces would result 


for any other shaped cavity. 
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However, only partialt consideration has been given to those 
points which lie in the coordinate planes separating the 
first octant from the rest of p-space. Obviously, a 
Significant number of points—those with zeros in their mode 
numbers—will be miscounted if Eq.(5.6) is not modified. 
Since these points represent the "unambiguous" mode numbers 
Cite, those with multiplicities less than two), a 
bifurcation of Eq.(5.6) is necessary in order to properly 


eoune —DOtENeITM and TE modes: 


By considering only the first octant of p-space, the 
off-axes points on the coordinate planes have effectively 
been cut into halves, and the points on the coordinate axes 
have DeenmeCcULS® INntOQesLourths. Followingmethis s:linemror 
reasoning, the modified forms of Eq.(5.6) are obtained as 


follows: 


The first quadrants of the ellipses formed by the 
intersection of the coordinate planes and the ellipsoid 


(Se2ehave areas 
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tIn this volumetric representation of characteristic points, 
it is understood that a point has not been fully counted 
until an entire equivalent unit volume, which can _ be 
associated with it, has been accounted for. It TS9Sinewwe this 
sense that "portions" of a lattice point can be counted 


separately. 
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and the line segments formed by the intersection of the 


coordinate axes and the ellipsoid (5.2) have lengths 


Usteh) 


in the first octant. The terms (5.7,5.8) can be converted 
to volume terms by multiplying by a simple numerical 
fraction representing the portion of volume required. Now 
the set of TE modes includes those represented by 
characteristic points in the p,p3; and p2p; planes, but 
excludes those represented by characteristic points in the 
P:iPpz2 plane; and conversely for the set of TM modes. The 
characteristic points on the coordinate axes represent 
longitudinal modes which are excluded from both TM and TE 


sets. Thus 


N+4,(w) = Now) = (S ayes 3) w? a Gist BomLit) Ww (5.9) 
167 Ca 4n c 

N+_(w) = No (w) te (Spams) w? = (L,+L2+3L;3) Ww G5r,.1'03) 
167 Ce an C 


RHeLeES wee eb b4+2loL, andseSs—2byl7.8m ties plUS7 MinUSm—S1GUS 
on the RH sides are required in order to correct No(w) for 
the undercounting/overcounting of the characteristic points 


described above. 


One further modification term in each of Eqs.(5.9,5.10) 
is required in order to properly account for the 000 mode at 
the origin of p-space. That is, each of the terms on the RH 


sides of Eqs.(5.9,5.10) contribute a positive or negative 
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fraction of the 000 characteristic point to their respective 


"subtotals". The sums of these contributions in each case 


are 


coo lech We eae es. (5.11) 
8 8 8 8 8 8 8 8 
due to Nt,,(w), and 
ingelda t+ >t - 1 + ar 4 (5.12) 
8 8 8 8 8 8 8 8 


due to Nt- (w), and these must be subtracted from 
Eqs.(5.9,5.10) respectively. Hence the partial asymptotic 
Spectral distributions of TM and TE modes respectively ina 


rectangular cavity are 


Ny y4(w) = No (w) = US eeaoa) w? av (GatleeLoy) Ww =P ») C5 re i3)) 
167 Ce 4n e! 8 


Ny _(w) = No (w) uF (Sen) w? FF (GE hobe cresyl) Ww oP a (Sina) 
167 oe an c 8 


and therefore the total asymptotic distribution of 
electromagnetic modes is 


w) = 2No (ai) a (Leto tine) Ww me ul ays 1s)), 


Newat 
= 27 Cc Z 


which is equivalent to the form (in terms of k) derived by 
Ralresmeecteal.. (27,2917 = andsewetys) close to the form 


anticipated by Pathria [23] in 1966. 


From the partial distribution functions (5.13,5.14), it 


is clear that any calculation of N—either analytical or 
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numerical—that presumes TM and TE modes to have the same 


distributions, will be in error. 
5.3 Asymptotic Spectral Density of EM Modes 


Applying the same reasoning which led to Eq.(2.51), the 
number of electromagnetic modes with resonant frequencies 


between w and wtdw is 
Deylw) dw = Ney(wtdw) - Neyy(w) (5.16) 


which, in terms of the p-space model, corresponds to _ the 
number of characteristic points within an ellipsoidal shell 
which is nou in general, of uniform thickness. 
SUDStICUCINOgMmEOS 1(5)0, 5.010) mm INtOeBEG .( Dlojmandmpertorming 


the subtraction yields 


Depylw) dw 2) arate) Oa 
Sige fel’ 
Sith tov lel otow) 0 (Se 1) 
27 c 
or 
Deylw) dw = V_ 1 [3w?bwt3w(dw)?+(dw) *] 
Shige el 
a (Ligtbo tla) a bw (5G) 
20 ©. 


which, as Sw approaches an infinitesimal, reduces to 


Deyy(w) dw =o Doo) 0Om a (lint bisa Oco (5.19) 
27 Cc 
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where Do(w)=dNo (w) /dw. The asymptotic spectral density 


function for electromagnetic modes in rectangular cavities 


is therefore 


Deylw) = 2Do(w) - LysU th) 4 (5220) 
Tv 


c 


which, as anticipated in Chapter 2, is just the derivative 


GEE. Gogelomwithe respect «too. 


The asymptotic forms for the component partial spectral 
densities can be derived in exactly the same manner as 


Eq.(5.20) above. The resulting formulae are 


Dryloum= Dolo)cC SiS 7 -Ss) OUE-BN LY tLe tb) Gt CSe2n) 
81 ce 4n © 

Dp (w) = Da (o) a7 (ins) w ee (L,+L2+3L;) au (53°22)) 
8n ce 4n (e 


which are the derivatives of Eqs.(5.13,5.14) respectively. 


5.4 Verification of the Asymptotic Formulae 


Onlvyethe wiirSt term inseach Of thesasymptocicmetormulee 
Survives in the limit wV7~, which is to be expected if 
Eq.(5.20) is to be consistent with the familiar black body 
radiation laws. These laws are well-tested experimentally 
and as such the asymptotic formulae are verified in the high 
frequency limit. However, at lower frequencies when the 
contributions of the second (and third) terms become 


important, Eqss(5. 1595.20) have yet to be verified 
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experimentally. 


The accuracy of the asymptotic formulae at lower 
frequencies was investigated by comparing the numerically 
and experimentally determined mode distributions of Chapter 
4 with the corresponding results predicted by the asymptotic 
formulae. This comparison is illustrated in Figure 5.1. for 
Ehesespectral distribution) of modes, and) in Figure 5e2 for 
the spectral density of modes. Also plotted on Figures 
SAI ea’ are the incorrect distributions predicted by 


Eq.(5.5) and its derivative respectively. 


A visual inspection of the curves in Figures 5.1,5.2 is 
quite convincing, but in order to provide some quantitative 
measure of the accuracy of the asymptotic approximations, an 
average error is quoted for each of the asymptotic curves. 
The figure E is defined for either N or D type (N,D) 


distributions by 


‘ 
Enrp = 4 CNS Deis - N,D(w, ) dw (Gays #2553) 


i=1 


where P is the total number of points calculated. The 
average error is useful for comparing the "fit" of different 
approximations, but does not give any indication of the 
amplitude of fluctations in the exact curve, about the 
approximate curve. A Slightly different figure—the average 


absolute error—is defined for either N or D type 
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Figure 5.1 Comparison of the exact and approximate spectral 
Gistribution of EM modes for a particular rectangular 
Cavity. The interior dimensions of the cavity are L,=23 cm, 


L2=50 cm, L3=52 cm. The dotted curve is the erroneous 


Gistribution, predicted) by Eqn (5.5). 
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Figure 5.2 Comparison of the exact and approximate spectral 
density of EM modes for a particular rectangular cavity for 
a bandwidth of 50 MHz. The cavity dimensions are the _ same 
SCmeLOtm rigiure tS. 1. The dotted ‘curve 1s the erroneous 


distribution predicted by the derivative of Eq.(5.5). 


5 7 


aisle | 


73 


distributions by 


PB 
Anrp _ > LCN *;D) lc om N,D(w; )6w| (5.24) 
=1 


and can be used to provide this extra information. 


The singular case presented in Figures 5.1,5.2 does not 
clearly show the effects of cavity shape and symmetry on the 
accuracy of the asymptotic formulae. In order to get an 
indication of these effects, Eqs.(5.5,5.19) were used to 
generate the approximate mode distributions for each of the 
hypothetical cavities discussed in Chapter 3. These results 
were not graphed, but the figures E,A for each case are 


compared in Table 5.1. 
5.5 Exact Spectral Density of EM Modes 


BCS oe 20-52 2) are smooth, continuous, monotonic 
functions of frequency and as such are only approximations 
for finite (if not too small) frequencies. However, if the 
bandwidth and/or the cavity is very small, then the exact 
Gistribution Dt can exhibit large fluctuations about the 
asymptote D(w) due to the step discontinuities in the exact 
distribution Nt. Of course the accuracy of Eqs.(5.20-5.22) 
is seriously degraded under these conditions. The exact 
spectral density of modes is, however, known analytically 


for rectangular domains and is given in [31] as 
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Table 5.1 Comparison of the figures E and A for the asymptotic 


formulae applied to the hypothetical cavities discussed in 
Chapter 3. The first figure of each pair is E as defined by 


Eq. (5.23) and the second figure is A as defined by Eq. (5.24). 


E AND A FOR CURVES OF D vs w/w, 


OF FOR 
APPROXIMATE | CURVES OF ; 
FORMULA N vs w/w, | BW=50 MHz | BW=100 MHz | L,=L,+lcem 
CUBIC CAVITY (L,=L,=L,=L=25cm) 
correct EM? | 0.00;2.19 | -0.14:2.70 | -0.12;3.14 | -0.03;2.67 
Necusnunerye |PEYAINZY al || blyaarce our | ocean) | ochikagsl sas 
RAYLEIGHx2° |-4.30;4.43 | -0.27;2.74 | -0.37:3.46 | -0.29;2.67 
RECTANGULAR CAVITY (L,=25cem, L,=35cm, L,=45cm) 
correct Em? | 0.02:1.63 | -0.01;2.08 | -0.03;2.50 | -0.04;2.95 
MCOUSTICNIGNIE 5586755260 | 2860; 3818 ule=5e25 Seco mine Seer once 
RAYLEIGHx2? |-5.68;5.68 | -0.19;2.09 | -0.38;2.53 | -0.39;2.98 


a I ne FF 


*mhe bandwidth associated with this column of figures is 100 MHz. 


these values were obtained using data generated from the correct 
asymptotic formulae for electromagnetic modes, Eqs. (5.15,5.19). 


“These values were obtained using data generated from the 
incorrect formula, Eq. (5.5), and its derivative. 


d these values were obtained using data generated from just the 
first terms of Eqs. (5.15,5.19). 
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wherenfij =p, c/4bj,0p'-(,07+f 2L2+f£202* and) the primessonmtness 
signs indicate that the f£,=0 and m=0 terms are to be omitted 
from the sum. ie is oe oe =e er a cube shaped cavity), 
then Eq.(5.25) reduces [31,32] to 
by +o 
Dewikle= Vk? | 1 + sin(2£kL) 
2f£kL 


f£,,f2,f£3;=-@ 


+o 
= IS | et ; '' cos (2mLk) | ton] 6. Gk) C526) 
27 2 


mMm=-o 


where f?=£7+f£3+f3 and the primes have the same meaning as 


Note that Eqs.(5.25,5.26) describe distributions in the 
mathematical sense, rather than distribution functions 
per se, and as such are valid only in integral form. Thus, 
the convolution 

© 
| Dk) Wk Kg) ke (527) 


O 
yields the number of modes within the bandwidth described by 


the window function W. A thorough and rigorous treatment of 


the mode density fluctuations is given in [25]. 
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5.6 Discussion 


The utility, at least for rectangular cavity design 
purposes, of the asymptotic formulae is apparently 
questionable at best since they clearly do not provide 
realistic descriptions of the exact mode distributions. The 
Gonaw new Ol mera GUneS m5. eons TOTmECOULSCEE TA TEENEGLOMEE tie 
asymptotic regime, and hence one could only expect "best 
fit" representations of the exact distributions. This may 
be adequate for large cavities or large bandwidths where the 
computation of the exact N,D is not feasible. However, for 
a large class of microwave cavity problems, the lower 
frequency range is the most interesting [42-49] and accurate 
knowledge of the spectral density of modes is _ desirable, 
especially ne bandwidth is limited. The approximate 
distribution functions would seem to be of limited value in 
these applications. Yet the approximate formulae (albeit 
incorrect versions of them) have been mentioned several 
times in the engineering literature, without qualification, 
in the context of microwave cavity design. But quite 
unexpected results could be obtained if Eq.(5.15), or 
especially Eq.(5.20) were used naively in design 


calculations. 


For nonrectangular cavity design, however, the 
approximate formulae may be very useful. As mentioned in 


the last chapter, there may be advantages to using 
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odd-shaped nonrectangular cavities for some microwave 
heating applications. In such cases, it can be extremely 
difficult, if not impossible, to calculate the exact mode 
astern routs on; and hence the asymptotic formulae yield 
valuable information, even for low frequencies. 
Furthermore, the accuracy of the approximations would be 
better for the case of irregular cavity shapes because of 


the absence of degeneracies due to geometric symmetry. 
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CHAPTER 6. SUMMARY AND CONCLUSIONS 


The purpose of this chapter is to briefly summarize the 
contributions made by this work to the field of Electrical 


Engineering. 


6.1 Regarding Calculated Mode Distributions 


A numerical algorithm has been presented which 
correctly enumerates the TM and TE modes within a_ finite 
bandwidth for an empty rectangular cavity. Three computer 
programs which utilize this algorithm for the 
design/analysis of microwave heating cavities were also 
presented. These programs are useful for analysis purposes 
and also for obtaining target cavity dimensions (i.e. a 
Starting point for the design process) which yield 


favourable mode distributions for a particular application. 


One of the programs was used to generate theoretical 
(exact) mode distribution data for several hypothetical 
cavities. The results revealed certain characteristics—in 
particular the marked sensitivity to cavity shape and 
bandwidth—which suggest that, due to pracnical 
considerations, intrinsic mode density calculations are 
potentially misleading from a design point of view. 
However, these calculations are still valuable for reference 


purposes, and additional numerical results were compiled in 
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reference tables which are intended to replace (and augment) 
previously published incorrect data [59]. The discrepancies 
between the correct and incorrect mode counts were shown to 
be on the order of 50-100% and hence it is extremely 
important to use the correct algorithm, especially when 


small cavities (or low frequencies) are being considered. 


COorrec teeasymptotic. wmode §distributvone eLunctvoOnSsmet oc 
rectangular cavities have been derived in a heuristic manner 
which can be understood by engineers without a 
graduate-level background in mathematics. These derivations 
emphasize the origin and nature of the physical differences 
between TM and TE modes which must be understood if these 
functions are to be used for research (or design) purposes. 
The formulae quoted here appear for the first time, in their 
correct forms, in the microwave engineering literature. 
Several factors affecting the utility of these analytical 
formulae were evident from comparisons with results obtained 


numerically and experimentally. These were discussed. 


Only rectangular cavities have been considered in 
detail here. However, it is understood that the physical 
principles which apply to mode counting in rectangular 
cavities also carry over to other geometries. For example, 
the fact that TM and TE modes have different distributions 
and therefore must be counted separately in order to obtain 


the correct number of modes in a cavity, holds true for any 
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regular geometry (i.e. where the distinction between TM and 
TE modes makes sense). Also, conclusions regarding factors 
which affect the mode distributions (e.g. shape, bandwidth) 
in rectangular cavities can probably be generalized, at 


least qualitatively, to other cavity shapes. 


6.2 Regarding Measured Mode Distributions 


The low frequency spectral distribution of 
electromagnetic modes in an empty, rectangular, laboratory 
cavity has been determined experimentally. These 
measurements constitute the first experimental verification 
of the electromagnetic mode counting theory [60]. The 
Straightforward meaSurement technique described here is 
applicable to microwave cavities of arbitrary shape and is 
therefore useful in investigations of nonrectangular 


Cavities and cavities with mode stirrers, turntables, etc. 


The effect of the presence of a mode stirrer on _ the 
Spectral distribution of modes in the laboratory cavity was 
investigated experimentally. The results indicate that a 
mode stirrer causes a more thorough splitting of degenerate 
modes, where the actual resonant frequencies depend strongly 
on the angular position (and, of course, the geometry and 
location) of the mode stirrer. Lower Qs and a very strong 
dependence of the coupling on angular position (etc.) were 


also observed. All sofe. theses, eitectse Cables Cc mUSCOmsCO 
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advantage in the design of microwave applicators. However, 
more theoretical and experimental work is required in order 


to exploit these effects optimally. 


6.3 Some Topics for Further Research 


In achieving the objectives set forth in Chapter 1, 
several interesting research topics have been uncovered 
which are unfortunately beyond the scope of this. thesis. 
For example: mode distribution and field pattern 
measurements on nonrectangular (e.g. trapezoidal) cavities; 
experimental (and theoretical) investigations of cavity 
perturbations due to mode stirrers; more thorough and 
rigorous investigations of the cross-coupling between 
multiple source antennae; corrections to N(w) to account for 
the presence of a mode stirrer; improvements to CAD programs 
for cavity design. These are some of the topics which will 
soon become very important in the area of small, low power 


applicator design (e.g. solid state microwave ovens). 


6.4 Concluding Remarks 


At this stage, only a qualitative interpretation can be 
given to the intrinsic mode distributions for a given 
cavity. Hence, it would seem that less (but not zero) 
emphasis should be placed on mode counting as a design aay 


at least for cavities intended for use aS microwave power 
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applicators. However, as research (or diagnostic?) tools, 
mode distribution measurements and calculations may be very 
valuable. Understanding these distributions is an essential 
first step in the investigation of shape and symmetry 


effects, and arbitrary loading in small cavites. 


Furthermore, as electrical engineers become more 
involved mn work which requires a knowledge of 
thermodynamics and radiation physics (e.g. as in the space 
program, fusion research, solid state electronics, etc.), it 
will be necessary to understand and apply the basic laws of 
Planck, Wien, Stefan-Boltzmann, Einstein, etc. The 
asymptotic spectral density formula is involved in all of 
these laws and the classical texts do not consider the 
corrections which become important at lower frequencies 
(ime smo intrared).. With a correct understanding of the 
theory of mode counting, it is less difficult to determine 
the limitations of the classical theory, and therefore avoid 
the kind of errors which currently exist in the microwave 


engineering literature. 


The last text on microwave power engineering has 
certainly not been written, and every new text will 
undoubtedly deal with the subject of mode counting in some 
form or another. It is hoped that this work will serve, not 
only to correct the current literature, but also to provide 


future authors and researchers with adequate reference 
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material to prevent further errors from OCCunmingE tees 
also hoped that this thesis will inspire further research in 


this area, especially on the topics given above. 
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APPENDIX A. COMPUTER PROGRAM LISTINGS 


A.1 NEWCAV 
1 SRUN *WATFIV T=10 P=99 5=*MSOURCE* 6=*MSINK* 7=-OUT(*L) PAR= e 
a pest seats T=10 NOWARN NOEXT NOLIST | ee 
4 C ET LT TO TE SS I ET AEN 
. : | Located in MTS file: NEWCAV Revised: 29 September 84 
7 ( SO ST SE rr eS 
8 C es ee ore re 
9 cE MAIN PROGRAM 
10 ce S| 
i) C 
2 (e 
13 G AA - base dimension along the X axis (cm) 
14 C AS - step size for the increments of dimension 4 (cm) 
15 C AN - number of steps taken in the direction of A 
16 c BB - base dimension along the Y axis (cm) 
7 (C BS - step size for the increments of dimension B (cm) 
18 C BN - number of steps taken in the direction of B 
19 C DD - base dimension along the Z axis (cm) 
20 C DS - step size for the increments of dimension D (cm) 
AS e DN - number of steps taken in the direction of D 
ee C F1 - minimum frequency in the band of interest (MHz) 
AS Cc F2 - maximum frequency in the band of interest (MHz) 
24 (€ RE - response to indicate whether or not to rerun 
P43) C LR - last chance to prevent termination of program 
26 (C SR - response to indicate whether or not to save results 
ALT C LS - last chance to save results of current run 
28 C RC - response to indicate whether or not to change data 
29 Cc 
30 C 
31 C 
32 € OD - indicates what output device is currently in effect 
33 G EF - indicates that input errors have been detected 
34 C FR - remains TRUE only until the end of the first run 
355 Gg BC - counts the number of unsuccessful runs 
36 C SW - switch to indicate whether or not to reset BC 
Si C UC - counts the number of elements in CAV 
38 C VC - counts the total number of TE and TM modes 
39 C HCP - contains TM mode totals for each element of CAV 
40 C ECP - contains TE mode totals for each element of CAV 
41 C PMX - contains the maximum values of the input parameters 
42 (C OKP - indicates whether or not input parameters are ok 
43 C CAV - contains dimensions of the most mode-dense cavities 
44 e MTE - contains indices of the TE modes supported by above 
45 C MTM - contains indices of the TM modes supported by above 
46 c FTE - contains resonant frequencies of each mode in MITE 
47 g FIM - contains resonant frequencies of each mode in MIM 
48 C PTE - characters required to construct TE mode spectrum 
49 C PTM - characters required to construct TM mode spectrum 
50 G LBL - contains mid-band frequencies for labelling spectra 
51 C 
ey C 
ei RELUSGEAAMAS BEY bo EU RUS nE lines 
54 + CAVG1 09s bse HO; One bal Niet 05.503), 
a) ce PMX U7 2H ECU LEQ ee Opt e4 hes 1b 27 
56 INTEGERS 2eES Ui Ky Dia Nak cite sion Carvin 
57 + ANT BN DNS BCUCAVG peG.HC, SWa Al Ey MES 
58 1 LBL(11),MTE(10,50,3),MTM(10,50,3),ECP(10),HCP( 10) 
59 INTEGER*4 OD 
60 LOGICAL*1 EF,FR/.TRUE./,OKP(6) 
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NEWCAV Cont'd. 


61 CHARACTER*1 PTE(10,10,71),PTM(10, 10,71) 

62 CHARACTER*4 RE,LR,SR,LS,RC 

63 COMMON /REAL/ AA,AS,BB,BS,DD,DS,F1.F2 

64 COMMON /INTGR/ AN,BN,DN 

65 ON ERROR GOTO 22 

66 C 

67 é a rn aoa eran es ee wer ge ee a ee 
68 G An introduction to the program is displayed at the start 
69 (e of every session. 

70 C 

v4 G 

he WRITE(6,611) 

gfe) (c, 

74 C 

ee fe BC is set to 0 initially and following a successful run. 
hed C 

78 11 BC=0 

79 Swe 1 

80 G 

81 iC 

82 G Current values of AA,BB,DD are displayed (if not first 
83 C run! and user is asked whether changes are required and 
ee : if so is prompted for the new values. 

BE Cc 

87 12 LE SNOT AERP THENSDO 

88 WRIMEG = 2h) maa AA array 

89 WRITE(6,621) ‘BB’ ,BB 

90 WRITECTCE¥O2i0e DD" SDD 

91 WRITE(6,624) 

92 READS, O11 GeND=22, ERREZ 2 RE 

93 EC eG. IMO?) ORR Gcl2).” ie’ t) GIR RG ole@hs “IN! 7) 

94 + OR. URE JEG. ty? ) OR, RCO, 7 @" |) Ge) Ta) ae: 

95 END olf 

96 WRITE(6,612) PMX(1) 

97 READ(5,*,END=22,ERR=22) AA,BB,DD 

98 C 

99 G 

100 C Current values of AS,BS,DS are displayed (if not first 
101 C run) and user is asked whether changes are required and 
102 C if so is prompted for the new values. 

103 c 

104 G 

105 13 Dev NOt. FR )THEN LODO 

106 WRITEC 6,621) 2AS AS 

107 WRITE(6,621) ‘BS’ ,BS 

108 WREREGG. 62 19) “DS DS 

109 WRITE(6,624) 

110 READ(5,511,END=22,ERR=22) RC 


TECURCAEO RENO) VOR. (RCUEO> ono). OReGRC.EQR aN) 
+ BOR URCSEQs In JOKE RC; EO a Das) GOTO e14 
END IF 
WRITE(6,613) PMX(2) 
READ(5,*,END=22,ERR=22) AS,BS,DS 


Bn i ne el ek a rr 
Current values of AN,BN,DN are displayed (if not first 


run) and user is asked whether changes are required and 
if so is prompted for the new values. 
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NEWCAV Cont'd. 


121 
122 
123 
124 
125 
126 
127 
128 
129 
130 
T31 
132 
133 
134 
nso 
136 
rod 
138 
139 
140 
141 
142 
143 
144 
145 
146 
147 
148 
149 
150 
tot 
oa 
AOS 
154 
she) 
156 
View! 
158 
159 
160 
161 
162 
163 
164 
165 
166 
167 
168 
169 
170 
qe | 
dite 
Des 
174 
a 
176 
177 
178 
179 
180 
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14 IF(.NOT.FR)THEN DO 
WRITE(6,622) ‘AN’ ,AN 
WRITE(6,622) ’ BN’ ,BN 
WRITE(6,622) ‘DN’ ,DN 
WRITE(6,624) 


READ(5,511,END=22,ERR=22) RC 


PRO RC EQ” NO") .ORet RGsEO eno: 


WIORS HOR GRE Qa Nis) 


rf SOR. VREGEO RM nimi ORS RG&E Oe 20sGOr TOmi5 


END IF 
X1=PMX (3) 
WRITE(6,614) XI 


READ(5,*,END=22, 


ERR=22) AN,BN.DN 


Current values of F1,F2 are displayed (if not first run) 
and user is asked whether changes are required and if so 
is prompted for the new values. 


15) Baten. FR) THEN DO 
WRETEUGs623ney Eile hil 
WRITEUG., 623) GFZ @eri2 
WRITE(6,624) 


+ 


READ(5,511,END=22,ERR=22) RC 
LECCRES E OR NO OIEOR URC. COs Og Ones RC EG. Nas 


TORS CRC. EQns nM wcORa (ROG e Ose Os Gon logs 6 
END IF 
WRITE(6,615) PMX(4),PMX(5) 


READ(5,*,END=22,ERR=22) F1,F2 


Input parameters are passed to subroutine PARCHK which 
sets certain elements of OKP=.FALSE. if incompatibilities 
are detected. The elements of OKP are then checked and 
appropriate error messages are displayed corresponding to 
elements of OKP which are .FALSE. The user is allowed 2 
attempts at a successful run before execution terminates. 


a EE NEUEIE SEIS 


16 


17 


CALL PARCHK (PMX, OKP ) 


ER=aPAESE.. 
DOM PES 15.6 
PECENOTSOKP CE) )EFSc TRUE: 
CONTINUE 
LEV EG REN. DO 
BC=BC+1 
TPUBC.BO 73 HENS DO 
WRITE(6, 821) 
GOUTOESS 
END IF 
WRITE(6,811) 
IF(.NOT.OKP(1))WRITE(6,812) 
IF (.NOT.OKP(2) )WRITE(6,813) 
IF (.NOT.OKP(3) )WRITE(6,814) 
1F(.NOT.OKP(4) )WRITE(6,815) 
IF (.NOT.OKP(5) )WRITE(6,816) 
IF (.NOT.OKP(6) )WRITE(6,817) 
GOT 0F21 
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NEWCAV Cont'd. 


181 
182 
183 
184 
185 
186 
187 
188 
189 
190 
191 
192 
ese 
194 
195 
196 
197 
198 
199 
200 
201 
202 
203 
204 
205 
206 
207 
208 
2093 
210 
2t4 
212 
213 
214 
215 
216 
Pe 
218 
219 
220 
221 
222 
223 
224 
eco 
226 
227 
228 
229 
230 
231 
232 
ane 
234 
2a5 
236 
2c/ 
238 
239 
240 
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END IF 


18 


If input parameters pass the PARCHK tests they are passed 
to subroutine DISPRS where the allowable mode numbers 

and corresponding resonant frequencies are calculated for 
every possible set of dimensions. The sets of dimensions 
which describe the most mode-dense cavities are returned 
in the array CAV. The mode indices are returned in the 
arrays MTE,MTM and the corresponding resonant frequencies 
returned in the arrays FTE,FIM. 


CALE DISPRS(CAV\MTE,MIM,FTE,FIM,UC,HCP,ECP, VC) 
IF(UC.EQ.0)THEN DO 

WRITE(6,818) 

GO? T0120 
END IF 


Subroutine LINSPM generates a pair of normalized line 
spectra (i.e. one for TE and one for TM modes) for each 
cavity described by the elements of CAV. The normalized 
magnitudes are returned in the arrays PTE and PTM. The 
array LBL contains a list of evenly spaced mid-band 
frequencies for labelling the spectra. 


CALDBLINSPEMICAV]MTERMIM, FTE IMA UCORCRRECR Orie PEM, LBL) 


The printed output includes a listing of input parameters 
followed by a list of the supported modes and the two 
line spectra for each of the cavities described by the 
elements of CAV. A copy of the printed output will be 
appended to the bottom of the temporary file -OUT if 
desired. 


OD=6 
WRITE(OD) 7a17 AATEBROD. AS BS DS,AN, BNIDN Fa FZ 
WRITE(OD,712) 
DG 1S l= 15.UG 
HC=HCP (1) 
EC=ECGP (1) 
WRITE(OD,713) 
TEQHCENE.O)THEN DO 
WRITECOD. 714) (COVit dy ewas 5 3) (MEI hy Kesey Sess) 
IF(HC.GT.1)WRITE(OD,715 
((MTE(I,L,M) 
IF(HC.EQ.50)WRITE(OD,81 
WRITE(OD,716) HC,’ TE’ 
PRCEC. NE. OJ THEN DO 
WRITE(ODS 71D im CUM IMAwLEs MEd MES sr iM La Ded be EO) 
IF(EC.EQ.50)WRITE(OD,819) 
END IF 
WRITE(OD,716) EC,’ TM’ 
EUSESDG 
WRITE(OD,714) (CAV(I,J), 
IF(EC.GT.1)WRITE(OD,715) 
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NEWCAV Cont'd. 


24) 
242 
243 
244 
245 
246 
247 
248 
249 
250 
25) 
252 
255 
254 
255 
256 
Jon 
258 
250 
260 
261 
262 
263 
264 
265 
266 
267 
268 
269 
270 
271 
oe 
203 
274 
aii 
276 
rat 
278 
279 
280 
281 
282 
283 
284 
285 
286 
287 
288 
289 
290 
ao 
292 
2Js 
294 
230 
296 
257, 
298 
299 
300 


DG GiGi iG iG) 


DIQ'S'@' @'@ 


20 


21 


((MTM(I,L, 
IFCEC.EQ. 50) WRITE (OD, 
WRITE VOD A716) ECs TM: 
WRITENODE 7 16)) HOW TE 

ENDS 

IF(HC.NE.O)THEN DO 
WRITE(OD,717) 
WRITE VOD 71807 CCRTECD Ni Pls Petage , N= 1610) 
WRITE(OD,719) 
WRITE(OD,720) 
WRITE(OD,721) (LBL(Q),Q=1,11) 

ENDS Ue 

IF(EC.NE.O)THEN DO 
WRITE(OD, 722) 
WRITE(OD, 718). ((PTIM(IGR,S)Se107 1) R=15 10) 
WRITE (OD,719) 
WRITE(OD, 720) 
WRITE (OD,721) 

ENDED 

CONTINUE 

IF(UC.EQ.10)WRITE (OD, 820) 


Mat Sie rN dele Leo Ge 
) 


CEBLUV in Ve, 11) 


The user is asked whether or not a duplicate copy of the 
printed output is required and if so the copy is appended 
to the bottom of the temporary file -OUT. 


IF(OD.EQ.7)WRITE(6,619) 
TEXOD SEQ SE7IHEN DO 
WRITE(6,618) 
READ(5,511,END=22,ERR=22) SR 
LEC (SR. EQ YESMINOR. SRO. tyesen. ORS PORE. we) 
JOR. USREEQ 2 ya MORAUSR EOhs OK Si OR Mt SREOR OK) 
SOR SUSR REO tec HENMDU 
OD=7 
GOATON 18 
END IF 
1F((SR.NE.’ NO’ ).AND.(SR.NE.’no’ ). AND. (SR.NE.‘N’ ) 
SAND. USRENE 2! m@)sANDee SRO NE" 008) THEN) DO 
WRITE(6,617) SR 
READ(5,511,END=22,ERR=22) LS 
1 P(ChES® E'QR” NOE E ORFS FEO noe OR LS Ne) 
FORCES SEOne fe One USO, 0's Gla C20 
END IF 
WRITE(6,620) 
END IF 


I 


Multiple reruns may be requested. The number of reruns 
is limited only by the maximum allotted CPU time. 


i 


SW=2 

WRITE(6,616) 

READ(5,511,END=22,ERR=22) RE 

TRCIRE EO. NES SOR SURE, 2 Ves JEORSIREZ EQ. 19) 


UORMURE EQ ™ Yul VOR SURE EUR UK MO Re RES EOS” Ob 5) 
+.0R.(RE.EQ.’ 1’) )THEN DO 


FR=. FALSE’. 
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301 GO FiO Bb1 2h) aSW 

302 END IF 

303 IF((RE.NE.’NO’).AND.(RE.NE.’no’). AND. (RE.NE.’N’ ) 

304 +. AND.(RE.NE.’n’). AND. (RE.NE.’0’)) THEN DO 

305 WRITE(6,617) RE 

306 READ(5 7511 ,END=22,ERR=22) LR 

307 TECUER SEQS NOS) OR. CERIEO!’ no's) (OR. CER SEOs? NU 

308 + ORL RGEOS ne) OR WER EOS Oa GOL Te. 24 

309 GO TO 99 

310 ELSE DO 

Sit GO TO 99 

a2 END IF 

a3 22 WRITE(6,822) 

Orr @ 

SS C 

316 C FORMAT STATEMENTS 

< y) C 

318 G 500 - 599 formatted input 700 - 799 program output 
aie 600 - 699 instructions 800 - 899 error messages 
321 G 

322 511 FORMAT(A4) 

323 Git EORMAIN/S 7 /// S111 10K 2 10, 

324 + ‘This program helps the user to design ', 

325) * ‘a rectangular multimode’/’ ' ,5X,’ cavity with the ', 
326 + ‘highest possible mode density given a particular set ’, 
Sai + /' ',5X,’of input parameters. These parameters ', 

328 = NCONS iSt Of seagseteofes tamting ?/ey 45% = dimensions 
329 + ‘(AA.BB,DD) and corresponding step sizes (AS,BS,DS)’, 
330 + ‘and’ /' ‘',5X,' step numbers (AN,BN,DN); and a pair ‘, 
331 + ‘of frequencies (F1,F2) to define’ /' ‘,5X,‘the ', 

See + ‘bandwidth of interest. The program then calculates ', 
SoS + ‘all the modes’ /' ‘ ,5X,' supported by each cavity ' 

334 + ‘formed by varying the starting dimensions.’ 

325 + /' ',5X,'’ Line spectra are displayed for the cavities’, 
336 - ‘ which yield the highest’/’ ' ,5X,’number of modes.’ ) 
NEY) 612 FORMAT(//'0' ,16X,’ ENTER BASE(STARTING) DIMENSIONS IN‘, 

338 me “CENTIMETERS 8) ees 2 4xie ADONGRX SY. ZOAXES RESPECTIVELY? -/ 
339 + ‘' ' 3X,’ [ the allowable range (of real numbers) ' 

340 + cise 10 SAA OEB.. BD) ge R yom cman’ 4) 

341 SUS ORM CH 1? ©" evs, TERNS Site Size OND (ELEN ME SIRS: UNE OK a4 
342 + SO TREGT IONS RESPECT VERY (fee s4Xwiitthe tal lowable ~ , 
343 + “range (of; realanumbers) 15:970 <sASsBS;DSecy 

344 + FJee ee ACU | me ieye 

345 614 BQRMAT (//"0': 912% SENTDERENUMBER OF YSTEPS. INeXPY¥ 2 DERECTIONS *, 
346 + ’RESPECTIVELY’/’ ‘' ,8X,’[ the allowable range (of ’, 
347 + binteoens | iSameOuc AN, BN DNC @oraec se 77) 

348 615 FORMAT(//’ 0’ ,8X,’ DEFINE BANDWIDTH BY ENTERING START,STOP ' , 
349 + “FREQUENCIES RESPECTIVELY’ /’ ‘,34X,’ IN MEGAHERTZ’ / 
350 + ‘ ' 4x,’ [ the allowable range (of real numbers) ', 

351 + 1p fete wee Ce Fe Gul GE ORD) alel Mi ae Be 

352 + ‘and F2-F1 < *°,F8.2,' MHz J*//) 

353 616 FORMAT(//‘0' ,26X,‘'DO YOU REQUIRE ANOTHER RUN?’ //) 

354 617 FORMAT(//' 0’ ,26X,’OK TO ASSUME ' ,A4,’ MEANS NO?’ //) 

B55 BiB SEORMAW//70. 9X." DOS VOU RWISH IOs SAVERIHE RESULT Se Crm ital ome, 
356 + “RUN IN THE FILE =OUT?"/7) 

357 619 FORMAT(//’0' ,33X,’ RESULTS SAVED’ ) 

358 696 FORMAT(//( 0’ ,ciX, RESULTS DISCARDED | 

359 624 FORMAN Ou, 21 eo7 HESCURRENT gVALUE DR emAg ey LSeabS. 2, cms 


360 622 FORMAT (2.0%, 2 ides THESCURRENTAVALUESOR Gen A2, teloe, Lo) 
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361 
362 
363 
364 
365 
366 
367 
368 
369 
370 
Sel 
372 
373 
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375 
376 
Sal 
378 
379 
380 
381 
382 
383 
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391 
392 
393 
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395 
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397 
398 
399 
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420 


DQ QO © G:@ 


623 FORMAT(’0’ ,21X,’ THE CURRENT VALUE OF ’ ,A2,’ IS’,F 12.4,’ MHz’ ) 
624 FORMAT(’ -' ,23X,'DO YOU WISH TO MAKE ANY CHANGES?’ //) 
711 FORMAT(’ 1’ ,20X,’ THE INPUT PARAMETERS USED IN THIS RUN ARE:’ 


+ PROS 2X, Air! Se Berea mies 0 een peetrer homme 
+ pmo eXo 4 Dl pzenES..28 im Cok mee Boker ASeria EB) 202 em: 
+ ings Bea Eos lia Cl? /upeeae tg w0Se eke. 2,0 cml 
+ pee 2a | ANgzi tal Suv bg 132 Ne Niwel el oye 22h Naa 
+ (0 ON A wes narrate Ty 
+ Wee LE op cieins (22 Caw Meza 
712 FORMAT('-',10X,'THE MOST MODE-DENSE CAVITIES ARE:’ | 
713 FORMAT(’ -' ,5X,64('_")/'0' ,8X,/ (cm)’ ,6X,/B (cm)’ ,6X,’D (cm)’, 
+ Baas hs (Mame 05,7 aE giMHZy 0 ee Xe Galle 1) 


TIS SRORMAT Sie, OA ol Rea AklienS 1 oe 2A bebo 
TASURORMA Lie Our ad Xe SLs Okie le a.) 


716 FORMAT (’0' ,5X,64("/ 3770’ ,6X,‘ TOTAL: ',13,1X,A2, 1X," MODE(S)’ 
+ TE Oe, 5X, BG lsartals 

717 FORMAT(‘-',10X,’NORMALIZED ELECTRIC FIELD LINE SPECTRUM FOR ', 
+ ETE MODES¢@7 7 Oe on 

718 FORMAT(’ *,5X,71A1) 

719 FORMAT(’ *.5xX,° ',10(' it 

ZO OMEORME Tiles 51a | eral Un camera |) a 

721 EORMAII. 209% 1115. 2kdy/0 0 see 

722 FORMAT(’-’,10X,’NORMALIZED ELECTRIC FIELD LINE SPECTRUM FOR ', 
+ PIM MODES: //20'5* 1) 

811 FORMAT('0’,15X,/*«*** ENCOUNTERED IMPROPER OR INVALID ’, 
+ ‘DATA amex’ /' ' 10X,!eex=ex MAKE CORRECTIONS TO ', 
+ ‘DATA FILE BEFORE RERUNNING ===»! | 

812 FORMAT(’0',10X,'’***«* ONE OR MORE "STARTING DIMENSIONS" ', 
+ ‘QUT OF RANGE «===! ) 


SiserORMAT Ue On ig see een ONes ORE MORES oer eol2o es ENDR OES ear 
+ ‘OUT OF RANGE *****' ) 
814 FORMAT(‘'0' ,8X,'’***** ONE OR MORE “NUMBER OF STEPS" ’, 
* ‘ENTRIES OUT OF RANGE ***»' ) 
815 FORMAT(' 0’ ,8X,’***** ONE OR MORE "START,STOP FREQUENCIES" ‘, 
+ ‘OUT OF RANGE **==! ) 
816 FORMAT(‘0’ ,10X,‘*=**** TOO MANY COMPUTATIONS REQUIRED TO ’, 
bs ~COMPLETE RUN *#****")) 
817 FORMAT(‘0O’ ,’***** DIMENSIONS TOO SMALL TO’, 


+ ‘SUSTAIN OSCILLATION WITHIN GIVEN BANDWIDTH **#**' ) 
818 FORMAT(' 0’ ,10X,’***** NO MODES WERE FOUND TO EXIST WITHIN ', 
+ ‘THE BANDWIDTH GIVEN *****' ) 
819 FORMAT(' 0’ ,10X,‘***** TOO MANY MODES EXIST WITHIN THE ', 
a: ‘BANDWIDTH GIVEN ****! /' / ,23X,'*##** OUTPUT LIST 
*: ‘TRUNCATED *****' } 
820 FORMAT('0O’ ,13X,’***** TOO MANY CAVITIES ARE EQUALLY ', 
Se PMODE=DENSE Sees! (03525) eee ee SOUL PUT Ll Sia, 
+ PTRUNCATED **¥#8%" ) 
B21 FORMAT(//'0' ,10X,’***** PLEASE FAMILIARIZE YOURSELF WITH ', 
+ ‘INPUT REQUIREMENTS ****"' /’ © ,21X,' "ee" BEFORE *, 
+ “ATTEMPTING ANOTHER RUN *****' //) 
822 FORMAT(//’ 0’ ,2X,’ ***** UNRECOVERABLE ERROR WILL TERMINATE ', 
+ “THIS RUN - PLEASE TRY AGAIN *****' //) 
99 STOP 
END 


Ss EUnet 


SUBROUTINE PARCHK 


Subroutine PARCHK provides some rudimentary error 
trapping by performing a series of simple tests on al] 
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421 C of the parameters (en mass) before allowing them to be 
422 C used in any computation. 

423 C 

424 C 

425 ¢ 

426 C AAM - maximum allowable value of parameter AA (cm) 

427 C ASM - maximum allowable value of parameter AS (cm) 

428 C ANM - maximum allowable value of parameter AN 

4293 C BBM - maximum allowable value of parameter BB (cm! 

430 6 BSM - maximum allowable value of parameter BS (cm) 

431 C BNM - maximum allowable value of parameter BN 

432 c DDM - maximum allowable value of parameter DD (cm) 

433 € DSM - maximum allowable value of parameter DS (cm) 

434 C DNM - maximum allowable value of parameter DN 

435 C F2M - maximum allowable value of parameter F2 (MHz) 

436 € BWM - maximum allowable value of bandwidth (MHz) 

437 C NCM - maximum allowable number of cavities 

hee : MNM - maximum allowable value of any L,M,N 

440 C 

444 C 

442 C OK - indicates whether or not input parameters are ok 
443 e AL - becomes TRUE if AA is less than one half wavelength 
444 C BL - becomes TRUE if BB is less than one half wavelength 
445 C DL - becomes TRUE if DD is less than one half wavelength 
446 Ec OKP - elements become FALSE if corresponding tests failed 
447 C 

448 C 

449 SUBROUTINE PARCHK (PMX, OKP) 

450 REAL*4 AAM,ASM,BBM,BSM,DDM,DSM,F2M, BWM, PMX(7) , 

451 + ALY ASY BB, ESOD, DSakiek2 

452 INTEGER*2 ANM,BNM,DNM,NCM,MNM,AN,BN,DN, I 

453 LOGICAL*1 AL,BL,DL,OKP(6) 

454 COMMON /REAL/ AA,AS,BB,BS,DD,DS,F1,F2 

455 COMMON /INTGR/ AN,BN,DN 

456 G 

457 C 

458 C All elements of OKP are initialized to .TRUE. and local 
459 G maxima are set according to values passed in array PMX. 
460 C sta Seon anti ett a i BE a a A ie a rt on 
461 c 

462 Dome Vestin6 

463 Th OKP UT) = 2 TRUE. 

464 AAM=BBM=DDM=PMX( 1) 

465 ASM=BSM=DSM=PMX/ 2} 

466 ANM=BNM=DNM=PMX (3) 

467 FQ2M=PMX( 4) 

468 BWM=PMX (5) 

469 NCM=PMX (6) 

470 MNM=PMX (7) 

471 G os ee ee eee ee 
472 C | 

473 C All input parameters must be equal to or greater than 

474 8 zero. 

475 C O_O 
476 C 

477 TELUAA LE. O)SCORTUDBSLE; GJe ORs LPO abies Os VORP inl) =a AL obs 

478 TE ( (AS@ETOO) ORS (BSH 0)" ORFUDSe LI. OF SOKR (2) eR FALSE. 

479 IF((AN.LT.0).OR.(BN.LT.0).OR.(DN.LT.0))OKP(3)=. FALSE. 

480 TEC (ETSLESO)F ORE (F20LE .0)MOKP 14) =. RALSES 
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481 
482 
483 
484 
485 
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487 
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494 
492 
493 
494 
495 
496 
497 
498 
499 
500 
501 
502 
503 
504 
505 
506 
507 
508 
5093 
510 
511 
5) 12 
Bis 
514 
31'S 
516 
517 
518 
519 
520 
ben 
22 
2S 
524 
525 
526 
527 
528 
529 
530 
531 
532 
53s 
534 
535 
536 
537 
538 
539 
540 


(PUMP NAMIE) GG) IGG) iIG) iG) (QAM) 


Oi GiGi GiGi) GCG GQ ©) OG GQ G@ 


All input parameters must be equal to or less than their 
respective maximum allowable values and must not combine 
in such a way as to cause an error. 


IF ((AA.GT.AAM) .OR.(BB.GT.BBM) .OR.(DD.GT.DDM) )OKP(1)=.FALSE. 
IF((AS.GT.ASM).OR.(BS.GT.BSM) .OR.(DS.GT.DSM) )OKP(2)=.FALSE. 
IF ((AN.GT.ANM) .OR.(BN.GT.BNM) .OR.(DN.GT.DNM) )OKP(3)=.FALSE. 
TEULE2.GT R2MJ. OR (GF 1 Gl. F2) ORs FO-F 1G 1. BWMoehOKPd 4)'= oF ALSE 
IF ((AN*BN.GT.NCM) .OR. (BN*DN.GT.NCM) .OR. (DN*AN.GT.NCM) .OR. 

+ (AN*BN*DN.GT.NCM) )OKP(5)=.FALSE. 
IF ( ( (AA+(AS*AN) )*2*F2*1E6/3E10.GT.MNM) }OKP(5)=.FALSE. 
IF (((BB+(BS*BN) )*2*F2*1E6/3E10.GT.MNM) )OKP(5)=.FALSE. 
IF (((DD+(DS*DN) )*2*F2*1E6/3E10.GT.MNM) )OKP(5)=.FALSE. 


The initial dimensions must be large enough to sustain at 
least one mode. 


PEVABSEP 2). Ute. 1E=107- THEN .DO 
OKPNG6) =. FALSE. 


GOTO. 9S 
END IF 
AL=BL=DL=.FALSE. 
TEVRATUAS*AN)™ LT ok 10/ (2b oe 1EO Alea TRUE. 
IF (BB+(BS*BN).LT.3E10/(2*F2*1E6))BL=.TRUE. 
IF (DD+(DS*DN).LT.3E10/({(2*F2*1E6) )DL=. TRUE. 
IF((AL.AND.BL).OR.(BL.AND.DL).OR.(DL.AND.AL))OKP(6)=.FALSE. 


If any elements of OKP become .FALSE. the user wil] be 
instructed (by MAIN PROGRAM) to examine the data file for 
improper entries. 


99 RETURN 
END 


SUBROUTINE DISPRS 


Subroutine DISPRS searches for those modes with resonant 
frequencies within the specified bandwidth - all possible 
combinations of 1,m,n within a calculated range are 
substituted into the dispersion relation for al] possible 
sets of dimensions. 


| See MAIN PROGRAM for a complete list of input parameters 


en nn nnn aU ISU EEE SEES 
i nn eae Utd dIE EEE ESNESESS SSS 


UC - counts the number of elements in the array CAV 

VC - counts the total number of valid TE and TM modes 

AC - current value of dimension ‘A’ used in computations 
BC - current value of dimension 'B’ used in computations 
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541 
542 
543 
544 
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548 
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550 
551 
sey 
SDs 
554 
Soo 
55S 
55y, 
558 
sie) 
560 
561 
562 
563 
564 
siete 
566 
567 
568 
569 
570 
ol 
eee 
oe 
574 
SfAe: 
576 
Shah 
578 
57 
580 
58 1 
582 
583 
584 
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586 
58 
588 
589 
S90 
Bo 
592 
os) 
594 
ail 
hele) 
eke 
598 
39g 
600 


DQ DWQOOQ@Q' QQ -e® 


ANININND 


ap Wate J ae oe IK pe 


GG G)iGiG) GCG) 


DC 

RF 

CAV 
MTE 
MTM 
FTE 
FIM 
(ome 
TT™™ 


resonant 
contains 
contains 
contains 
contains 
contains 
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current value of dimension ‘'D’ used in comoutations 


frequency calculated from dispersion rel. 
dimensions of most mode-dense cavities 
indices of the TE modes supported by above 
indices of the TM modes supported by above 
resonant frequencies of each TE mode in CAV 
resonant frequencies of each TM mode in CAV 


stores current values of L,M,N,RF for valid TE modes 
stores current values of L,M,N,RF for valid TM modes 


SUBROUTINE DISPRS(CAV,MTE,MTM,FTE,FTM,UC,HCP,ECP, VC) 

REAL SS SAAC AS, BEADS, OO ,US yr mre acs BOT UG, REV Aly aes 

Ie (SO 475 1 LN o0e 4) EITEUIO 50) FIM. 10, 50m CAV UO. | 
[INTEGER = 22 teu ks, CaM Noe OR AAN, BN ON SAS, AE De: Oba REs PER One 
UG CgeO VN CaVEn Vi nv DA Estat ty 

MIEIIO S05 Ss), MIMi TC SO; erm eCr (0) HCE yo) 

COMMON /REAL/ AA,AS,BB,BS,DD,DS,F1,F2 

COMMON /INTGR/ AN,BN,DN 


+ 


+ 
+ 


UC= 
VC= 


0 
0 


MAXIND(A1,A2,A3)=((A1+(A2*A3) )*2*F2*1E6/3E10)+1 


AE= 
BE= 
DE= 


AN+1 
BN+1 
DN+1 


PE=MAXIND(AA,AS,AN) 
QE=MAXIND‘BB,BS,BN) 
RE=MAXIND(DD,DS,DN) 


AC,BC,DC comprise the current set of dimensions. | 


DO 


coe T2417, Ae 


AC=AA+( (1-1) *AS) 


DO 


19 J=1,BE 


BC=BB+((U-1)*BS) 


DO 


1S) (Xe 4} 5 (DIE 


DCS B DGC i * 0S 1 


L,M,N comprise the current set of mode indices. 


VH=0 


VE= 
DO 
L=P 
DO 
M=Q 
0 
N=R 


Sy 


0 
1 dpe ak gk 
aa 
11 Q=1,QE 
ah 
i= RE 
vil 


AC,BC,DC,L,M,N are all substituted into the dispersion 
relation and if the resulting resonant frequency RF is 
within the specified bandwidth then L,M,N,RF are either 
stored (temporarily) in TTE,TTM or rejected depending on 
whether the mode is TE or TM or neither. 
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601 C 

602 RF=SQRT(((L/AC)**2)+((M/BC)**2)+((N/DC)**2) )*(3E10/2) 
603 PECRP SL pe (EI*4E6) GO) TOO 14 

604 LPURF. Gla tpe*1£6))GO TOs 14 

605 IF(((L.GT.0).OR.(M.GT.0)).AND.(N.GT.0))THEN DO 
606 VH=VH+ 1 

607 TEV Te 

608 TTE(VH,2)=M 

609 TTE(VH,3)=N 

610 TTE(VH,4)=RF/1E6 

614 TEA Vitee OO. SD) GOnTOmd2 

612 END IF 

613 IF((L.GT.0).AND.(M.GT.0)) THEN DO 
614 VE=VE+1 

615 TMV = 

616 TIM(VE,2)=M 

(5391 7 TTM(VE,3)=N 

618 TIM(VE,4)=RF/1E6 

619 LEM VETEO 350) GOMTORt2 

620 END IF 

621 11 CONTINUE 

622 C 

623 C 

624 € Dimensions of the most mode-dense cavities are saved in 
625 G the array CAV and the corresponding mode numbers and 
626 C resonant frequencies are saved in the arrays MTE,MTM and 
627 C FTE,FTM respectively. 

628 C 

629 € 

630 12 VT=VH+VE 

631 TEV Te BOR OiIGO Ome 

632 IF (VT-VC)19,14,13 

633 13 UC=0 

634 14) WCSUCse 7 

bio VC=VT 

636 HC=VH 

637 EC=VE 

638 HCP (UC) =HC 

639 IF(HC.EQ.0)THEN DO 

640 eet UC. 1 Se 

641 ns (We, I, HS 

642 Go) w@ Ws 

643 END IF 

644 ao) SS. ENS 7 lake 

645 PREAUG, AY Pe Le CH yea 

646 DOM Se N=atene 

647 15 MTE(UC,HY,HX)=TTE(HY,HX) 

648 iS j2C? (Ue see 

649 IF(EC.EQ.0)THEN DO 

650 BIM(UCy THe0 

651 MTM(UC,1,1)=0 

652 fa), y@) Alte! 

653 ENDE IF 

654 DOM t/a Gane 

655 FIM(UC,EY)=TIM(EY,4) 

656 DOM V/ae Ee N=minne 

657 17 MTM(UC,EY,EX)=TTM(EY,EX) 

658 ish (6 AAWHUE 5 WL) SENC. 

659 CAVCUC A 2 SBC 


660 CAV(UC,3)=DC 


ay 
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661 
662 
663 
664 
665 
666 
667 
668 
669 
670 
671 
672 
67s 
674 
675 
676 
677 
678 
679 
680 
681 
682 
683 
684 
685 
686 
687 
688 
689 
690 
691 
692 
693 
694 
695 
696 
697 
698 
699 
700 
701 
702 
703 
704 
705 
706 
707 
708 
709 
710 
(ag 
ne 
(Ae: 
714 
B15 
716 
Ti vit! 
718 
719 
720 
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GGG G)'G)' GG) 


105 


IF(UC.EQ.10)GO TO 99 
19 CONTINUE 


99 RETURN 


END 


SS 


SUBROUTINE LINSPM 


Subroutine LINSPM generates two (crude) normalized line 
specta for each cavity described by the elements of CAV. 


UC 

HCP 
ECP 
CAV 
MTE 
MTM 
ELE 
FIM 


wC 

XTE 
XTM 
WTE 
WIM 
TEU 
TMU 
TEN 
TMN 
Pais 
PTM 
LBL 


- counts the number of elements in the array CAV 


contains 
contains 
contains 
contains 
contains 
contains 
contains 


the TE mode counts for each cavity in CAV 
the TM mode counts for each cavity in CAV 
dimensions of the most mode-dense cavities 
indices of the TE modes supported by above 
indices of the TM modes supported by above 
resonant frequencies of each mode in MTE 
resonant frequencies of each mode in MIM 


MAIN PROGRAM for a complete list of input parameters 


counts number of modes with same X axis position 


contains 
contains 
contains 
contains 
contains 
contains 
contains 
contains 


the X axis positions of each TE mode 

the X axis positions of each TM mode 

the multiplicities of each TE mode 

the multiplicities of each TM mode 
unnormalized E field magnitudes of TE modes 
unnormalized E field magnitudes of TM-modes 
normalized E fieid magnitudes of TE modes 
normalized E field magnitudes of TM modes 


characters required to construct TE field spectrum 
characters required to construct TM field spectrum 
- contains mid-band frequencies for labelling spectra 


SUBROUTINE LINSPM(CAV ,MTE,MIM,FTE,FIM,UC, HCP, ECP, PTE, PTM;LBL ) 
REAER4RAASRS BEARS OD FDSHE U,.PeebrmnOs, COnCORHEREr, 


+ 


See oy 


CHARACTER* 1 


CAV‘ 10,3) 
DNGEGER*2 ei, Ii, 


FETE OFS GIS BIN (TORS OUTED (SOY SL MUES OH 
2, PORUEAUHAKE SR ANGN PO WPL eM eN, 


ANG BN ONG let ea Qe tree Moe Oromia 
HL,HM,HN,EL,EM,EN,UC,WC,EC,HC, 
MIEUTOSSU RS] MIME LONSO SIS AT EIS ORM oO, belie das 
WhE( 77) pWIMG7 7), ECP( 10) ,HCP (10); TEN(50),1 (50) 


PIE (105 10,74) -PIM(10) 10,219 


COMMON /REAL/ AA,AS,BB,BS,DD,DS,F1,F2 
COMMON /INTGR/ AN,BN,DN 


a 


The character arrays PTE and PIM are initialized with 
blanks throughout. The entire bandwidth defined by F1 
and F2 is resolved into 71 discrete X axis positions. 


en nn nnn SEU? 


vue oe 


. 
i é ; 7 
; an ee cAitni tere ati 7 
‘ ez +4 7) “i eT id “tts d > 
' as 7 Tt) - 
eed nel nt ee 7 o 


fh? gas. st ch /enies 

pa vi "<a > a Veg — | ‘ 
| Cah ae ytrvats dizer « Kite) See Ny oI ae = 
| gpbieves cirei~crn - ace Me ae | Ee é 


i | 
seria ed) tga -épor a! at) te a 
| game pe sakeis sian eh ar mi oe : ‘ 
ite + | @2 >be , Saree J pron art ae 
} GTB R46 WEem (ge Y Frau" an singe ur4q 


4 ee and 
, 


ed 
? 
: a? ea waa 
nein 
—— 4 a 
. — a : 
> J re | 
_ 
a an 
a 


ra 


| ned acd a 2b.8 mene eT im mea win’ So 7 7 
j - es ’ -_ ( 
ee = « 95 We’ eee, ore an 
tem 9% — ‘bts Bee : 
Ci 1 ciel yates 7 2 
: a -<-* 
7 yal 


4 
ae os id she a 
a eth 
| “er Nu?) © 
' me THe etay? wt a 
i 


my 7 
~—~— 


i 
- 

7 an, 

y a 


* 


Le, 77Sa5 Stas?) t 
ee a ja fieetee. 


7 =* , 4 ot 


. 
sahitedination o— rm See 
_ 


a? - 


GewUEecoveveviuumsugeuyavurvey 
a 
_ 


0) aTaes sua ae afi 


ae Fis 
an pees: fr 
? 
ie 


~ 


hrs 284 ig bare. ore 
-e a )weeieayy 


NEWCAV Cont'd, 


OO OG @iG) Gi 


PG? CRC WGI.C).€73 CG) CC CI) COC) 


iG 


12 


is 


14 


OUm ieleele levee 

DOS Agi 2 = 1810 

DO. 43 1351571 

PL Eat] 24 Loan ¢ 
CIMUDinieeiare 
CONTINUE 


The relative position along the X axis of each TE mode is 
located with respect to its resonant frequency. The 
number of modes found to occupy the same X axis position 
are counted by WC and stored in WTE. 


DOw 2 getie 1, UC 

HC=HCP (1) 

TE CHC.£0.01GO s10) 14 

DO 13 JH=1,HC 

ATE CUR = (CCRT EVI UR) =e 7 eo Sean ae Onin 
wC=0 

DO 12 KH=1.JH 

IF(XTE(KH) .EQ.XTE( UH) )WC=WC+ 1 
W=XTE( JH) 

WTE(W) =WC 

CONTINUE 


The magnitudes of the TE modes are calculated for each 
mode and stored (temporarily) in TEU. 


(lea 
(305 al 2 
(el eu Elns)} 
(Sleruitan) 


QBSCAVi 1,2! 

QD=CAV(I1,3) 
TEU(JUH)=TEFMAG(HL,HM,HF , QA, QB) 
CONTINUE 


The relative position along the X axis of each TM mode is 
located with respect to its resonant frequency. The 
number of modes found to occupy the same X axis position 
are counted by WC and stored in WIM. 


EC=ECP(1) 

PEREG EO] 07 GUST Dek 

DO 16 JE=1,EC 

XIMVE LS COLEIMGlSUE JQ GR iy (F2ckt)) #70) +4 
wC=0 

DOM IS ,KE=15 ce 

IF(XTM(KE).EQ.XTM(JE) )WC=WC+1 

W=XTM(JE } 

WTM(W) =WC 


15 CONTINUE 
| 
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be. C The magnitudes of the TM modes are calculated for each 
82 c mode and stored (temporarily) in TMU. TEU,TMU are then 
783 C passed to the subroutine MAGNRM which returns the 

oe : normalized magnitudes in the arrays TEN.TMN respectively. 
786 C Ee ee ek eT eer a ae geese eens 
787 EL=MTM(1,JE, 1) 

788 EM=MTM(1,JE,2) 

789 EN=MTM(1,JE,3) 

790 EF=FTM(1,JE) 

791 QA=CAV(1,1) 

792 QB=CAV(1,2) 

793 QD=CAV(1,3) 

794 TMU(JE)=TMFMAG(EL,EM,EN,EF,QA,QB,QD) 

795 16 CONTINUE 

796 17 IF(HC.NE.O)CALL MAGNRM(TEU,HC, TEN) 

797 IF(EC.NE.O0)CALL MAGNRM(TMU,EC, TMN) 

798 (€ 

793 c 

800 C The spectra PTE and PTM are formed by writing over the 
801 C appropriate blank characters with a plotting symbo! so 
B02 G that the resulting pattern resembles a plot of magnitude 
803 c versus frequency. 

804 C 

B05 C 

B06 TECHOSEOS OV GORTORAS 

807 DO 18 M=1,HC 

808 SH=TEN(M) 

B09 VH=XTE(M) 

810 PO) Wes Vi ieal, Sta 

811 Vee HOS yf ise | 

812 IO MATS WARE! (300) SUEMNE (aL NA MAR eS? 9h" 

813 TRGWILEGVH SC EOS2PPTECIAVE SVs" 2" 

814 IF(WTE(VH).EQ.3)PTE(1,YF,VH)=' 3’ 

815 IME MATE (MAa)) oleae qe (al Ae Mints? ct 

816 IM (Me (Mal) Qe its WC Nile Wnt) a? $5)" 

817 IE ATE (Wel) tel@). yeast, Ve Mini BY (Ss 

818 LECWLEAVEH) EQ. 7 )PTEUL YEU VR) =" 7" 

819 PRGWRECVH) EO... B) PTET YE, VA b= 8' 

820 IEE Ml) eek, SONS (ICG NAF 9 ni)? SP 

821 ISRIECWIEU VA) 2G) STPTECI YE Vins) = 

822 19 CONTINUE 

823 LPEEC. EC. 0) GUE TU sa | 

824 DOM ZOMNENEC 

825 SE=TMNIN) 

826 VE=XTM(N) 

827 DOPZORY2= lgiSE 

828 YP=10-Y2+1 

829 LECWIM Ves EOi 1) PLM eRe Pa Ga) = aad 

830 IF(WTM(VE).EQ.2)PTM(1,YP,VE)=' 

831 IF(WTIM(VE).EQ.3)PTM(1,YP,VE)=' 3’ 

832 IF(WTM(VE).EQ.4)PTM(1,YP,VE)=' 4' 

833 IF(WTM(VE).EQ.5)PTM(1,YP,VE)=' 

834 IF(WTM(VE).EQ.6)PTM(1,YP,VE)=' 6' 

835 IF(WTM(VE).EQ.7)PTM/1,YP,VE)=’ 7’ 

836 IF(WIM(VE).EQ.8)PTM(1,YP,VE)=’ 8’ 

837 IF(WIM(VE).EQ.9)PTM(1,YP,VE)=‘ 9’ 

838 20 IF(WTM(VE).GT.S)PTM(1,YP,VE)=' *’ 

839 21 CONTINUE 


840 DOM22) Pang Al 


20 8a NG ABN. - ae 


a 
a 
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22 
ee 


ag 


og 


LBL(P)=F1+((P-1)*(F2-F1)/10) 
RETURN 
END 


eS 


FUNCTION SUBPROGRAM TEFMAG 


Subprogram TEFMAG calculates the magnitude of the total 
electric field (to within an arbitrary constant) for a 
given TE mode. 


REAL FUNCTION TEFMAG(L,M,F,A,B) 
READS 4 RA BGKG, EAE Y 

INTEGER*2 L,M 
KC=((L/A)**2)+((M/B)**2) 
EX=(F*M)/(KC*B) 

EY=(F*L)/(KC*A) 

TEFMAG=SQRT( (EX**2)+(EY**2) ) 
RETURN 

END 


FUNCTION SUBPROGRAM TMFMAG 


Subprogram TMFMAG calculates the magnitude of the total 
electric field (to within an arbitrary constant!) for a 
given TM mode. 


REAL FUNCTION TMFMAG(L,M,N,F,A,B,D) 
REALH4 F,A,B,D,KC.EX,EY 

INTEGER*2 L,M,N 
KC=((L/A)**2)+( (M/B) **2) 
EX=(L/A)*(N/D)/KC 

EY=(M/B)*(N/D)/KC 
TMFMAG=SQRT((EX**2)+(EY**2)+1) 
RETURN 

END 


SUBROUTINE MAGNRM 


Subroutine MAGNRM normalizes an array YYU to the largest 
element in that array and copies the result to an integer 
array YYN whose elements have a maximum value of 10 and a 
minimum value of 1. 


SUBROUTINE MAGNRM(YYU,N,YYN) 
REAL*4 LV,YYU(50) 

INTEGER*2 1,U,N,YYN(50) 

LV=0 

DO 11 1 
DEY OH 
CONTINU 


=1,N 
1) 
E 
IF(LV.EQ.0 
( 
J 


.GT.LV)LV=YYU(1) 
LV=1 


(J)/LV)*10 
Q.0)YYN(J)=1 
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901 12 CONTINUE 
902 99 RETURN 
903 END 

904 /EXECUTE 

905 /BTCHEND 


ind of file 


pees: ru 


vines 


110 


A.2 MODIST 

. fk *PascalW T=25 P=50 SPRINT=-OUT(*1+1) PAR=NOLIST 

} 

3 { } 

: Located in MTS file: MODIST.NEW Revised: 09 August 84 

6 

if eae ModeDistribution(input,output); {$ S 1000000} , 

9 { } 
10 { Generates N or D vs. frequency data using the asymptotic mode |} 
11 {| distribution functions and compares the result with the exact |} 
He {| values obtained from the dispersion relation. 
14 

15 CONST pi=3.141592654; {computational value of pi} 

16 c=2.9979256E8; {speed of light in free space} 

i Tilt=1.0£-20; {smallest allowed denominator} 

19 VAR Lx,Ly,Lz: real; {cavity dimensions} 

20 fMin,fMax: real; {range of distribution} 

21 StepSize: integer; {resolution of distribution} 

22 bandwidth: real; {determines type of distribution} 

ce V,Sx,Sy,$z,L,R: real; {geometric properties of cavity} 

25 PROCEDURE ComputeVSL; 

26 { } 
27 } 
28 {| Assigns values to the (global) variables V,Sx,Sy,Sz,L. 
29 

30 

31 BEGIN {ComputeVSL} 

32 :=Lx*Ly*Lz; {volume} 

33 Sx:=2*Ly*Lz; {surface area of walls normal to x axis} 

34 Sy:=2*Lx*Lz; {surface area of walls normal to y axis} 

35 Sz:=2*Lx*Ly; {surface area of walls normal to z axis} 

36 L:=Lx+Ly+Lz; {sum of linear dimensions} 

37 END; {ComputeVSL} 

38 

39 PROCEDURE GenerateData; 

40 { } 
41 { } 
42 {| Generates exact and approximate N vs. f or D vs. f data } 
43 {| (and associated discrepancies) depending on the value of : 
44 {| bandwidth. 

45 { } 
46 ; 

47 VAR fOne,fTwo: real; {defines valid range for exact formula} 
48 fPrime: real; {moveable dummy abscissa} 

49 ErrorTE,Error™M,ErrorEM: real; feavacrs errors} 

50 ErrSumTE,ErrSum™M,ErrSumEM: real; {integrated errors} 

51 AbsErrSumTE,AbsErrSumTM, AbsErrSumEM: real; {int abs err} 
52 TE,TM,EM: integer; {exact number of modes} 

Se i,j,k: integer; {dummy indices} 

54 

55 PROCEDURE CountModes(fOne,fTwo: real); 

56 {I } 
ir 

58 {| Determines the number of allowable modes with resonant } 
59 {| frequencies between fOne and flwo. 
60 {Jo 
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61 

oe VAR px,py,pz: integer; {test mode numbers} 

ae eet MaxIndex(dimension, frequency: real): integer; 
66 } 
a Calculates the maximum value of a mode index. 
63 

i VAR HalfWaves: real; {dimension in terms of half waves} 
Ue BEGIN {MaxIndex} 

us HalfWaves :=2*dimension*frequency/c; 

74 MaxIndex:=trunc(HalfWaves ) 

ae END; {MaxIndex} 

a ASE TestMode(px,py,pz: integer); , 
79 { } 
80 {| Tests whether or not (px,py,pz) specifies an allowed |} 
81 mode and if so increments the appropriate counter(s). 
82 

83 

84 VAR RF: real; {resonant frequency} 

85 

86 BEGIN {TestMode} 

87 RF :=sart(sqr(px/Lx)+saqr(py/Ly)+saqr(pz/Lz))*c/2; 

88 IF ((RF>=fOne) AND(RF<=fTwo)) THEN 

89 

90 BEGIN {IF} 

g TFC (px>0)0R(oy>0) JAND(pz>0) THEN TE:=TE+1; 

92 IF (px>0)AND(py>0) THEN TM: =TM+1; 

93 EM:=TE+TM 

94 END; {1F) 

95 

96 END; {TestMode} 

97 

98 BEGIN {CountModes } 

99 FOR px:=0 70 MaxIndex(Lx,fTwo)+1 DO 

100 

101 BEGIN {FOR #1} 

102 FOR py:=0 TO MaxIndex(Ly,fTwo)+1 DO 

103 

104 BEGIN {FOR #2} 

105 FOR pz:=0 TO MaxIndex(Lz,fTwo)+1 DO 

106 

107 BEGIN {FOR #3} 

108 TestMode(px,py,pz) 

109 END; {FOR #3} 

110 

144 END; {FOR #2} 

142 

113 END; {FOR #1} 

114 

115 END; {CountModes} 

116 

WA FUNCTION NTE(frequency: real): real; 

118 ee} 
119 ; - 

120 {| Calculates NTE from the asymptotic (partial) formula. } 
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{| 


VAR Arg3,Arg2,Arg1: real; {reduced frequencies} 


BEGIN {NTE} 
Arg3:=exp(3*1n(frequency/c) ); 
Arg2:=sqr(frequency/c) ; 
Argi:=frequency/c; 
NTE: =(4*pi*V*Arg3/3)+(pi*(Sx+Sy-Sz)*Arg2/4) - 
((Lx+Ly+(3*Lz) )*Argi/2)i+( 3/8) = 
END; {NTE} 


FUNCTION DTE( frequency ,bandwidth: real): real; 


{ 
Calculates DTE from the asymptotic (partial) formula. 


VAR dArg3,dArg2, 
dArg!: real; {derivatives of reduced frequencies} 


BEGIN {DTE} 
dArg3:=3*saqr(frequency/c)/c; 
dArg2:=2*frequency/sar(c) ; 
dArgis=i/c: 
DTE:=((4*pi*V*dArg3/3)+(pix(Sx+Sy-S$z)*dArg2/4) - 
((Lx+Ly+(3*Lz))*dArg1/2))*bandwidth; 
ENO? {DLEt 


FUNCTION NTM( frequency: real): real; 
{| Calculates NTM from the asymptotic (partial) formula. 
{ 


VAR Arg3,Arg2,Arg1: real; {reduced frequency} 


BEGIN {NTM} 
Arg3:=exp(3*1In(frequency/c} ! ; 
Arg2:=sqr(frequency/c) ; 
Argi:=frequency/c; 
NTM:=(4*pi*V*Arg3/3)- (pix (Sx+Sy-Sz)*Arg2/4) - 
iia ais ewiouie pie 
END; {NTM} 


UNCTION DTM( frequency ,bandwidth: real): real; 


F 

{| Calculates DTM from the asymptotic (partial) formula. 

{ 

VAR dArg3,dArg2, 
dArg!: real; {derivatives of reduced frequencies} 


BEGIN {DTM} 
dArg3:=3*sqr(frequency/c)/c; 
dArg2:=2*frequency/sqr ic); 
dArg1:=1/c; . 
DTM: =( (4*pi*V*dArg3/3) -(pi*(Sx+Sy-Sz)*dArg2/4)- 
((Lx+Ly-Lz)*dArg1/2))*bandwidth; 
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END; {DTM} 


FUNCTION NEM( frequency: real): real; 
{ 


Calculates NEM from the asymptotic formula. 


VAR Arg3,Arg1: real; {reduced frequencies} 


BEGIN {NEM} 
Arg3:=exp(3*1n( frequency/c) ) ; 
Arg1:=frequency/c; 
NEM: =(8*pi*V*Arg3/3)-(L*Arg1)+(1/2) ; 
END; {NEM} 


ia DEM( frequency ,bandwidth: real): real; 


{| Calculates DEM from the asymptotic formula. 


VAR dArg3,dA4rg1: real; {derivative of reduced frequencies} 


BEGIN {DEM} 

dArg3:=3*sar(frequency/c)/c; 

dArgiczi/c; 

DEM: =((8*pi*V*dArg3/3)-(L*dArg1))*bandwidth; 
END; {DEM} 


PROCEDURE OutputNvsf: 


Tabulates the values of N computed from the approximate 
formulae and the deviations of each from the exact 
value of N. 


ae caetien tation hemetien Ramin amt) 
ot a ed ged et et 


PROCEDURE PrintTitles; 
{ 


} 
{ } 
{| Prints column headers for table. 
{ 


nn 


BEGIN {PrintTitles} 


Writeln (' N versus Frequency Data:’ ); 
Writeln; 
Writeln ( Ue 2 ; 

t ‘ ) : 
Writeln (' TH (EY el Cr TOY 
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Writeln ( 
Writeln ( 
Writeln ( 


‘ 
‘ 
‘ 
/ 
‘ 
‘ 
, 
/ 
‘ 
/ 
U 
‘ 


Writeln (' 


EXACT 
EXACT 
EXACT 


END: {PrintTitles 


flag PrintResults; 


i 9, 
‘ Q, 
ape Na 
APPROX 
APPROX 
APPROX 


. 
. 
. 
. 


ERROR ', 
ERROR ', 
ERROR |‘); 


i, laneter|meees ois. It): 
PET er ROME SD yey a Eee a 


{ 
{| Computes the total number of modes with resonant 
frequencies less than fTwo. 


BEGIN {PrintResults} 
Error le:=nE-NIE(fiwo): 
ErrorTM:=1TM-NTM(fTwo) ; 
ErrorEM: =EM-NEM(fTwo); 


hee U7 NV? 2 Payers WO), 7 7 We 
Mie Ue" le 3's." 
EGronitewes caw 
Wrate (Co [ee TM: 5. 
ErroriM: 7:2, 
Weviest oi, ENS 
ErrOGEMoa7 ten | a) 
Writeln; 


Enc Sumiles=erRsumMUETeEhhOoniie: 


ErrSumTM: =ErrSum?TM+Error tM; 
AbsErrSumTM: =AbsErrSumTM+abs(Error tM); 
ErrSumEM: =ErrSumEM+ErroreEM; 
AbsErrSumEM: =AbsErrSumEM+abs(ErroreEM) ; 
j:=j+1; {count the number of output lines} 
END; {PrintResults} 


BEGIN {OutputNvsf} 
IF (fPrime=fMin) THEN 


BEGIN {1F} 


PrintTitles; 


ErrSumTE:=0; 


ErrSumTM: =0; 
Err SumEM: =0; 
AbsErrSumTE:=0; 
AbsErrSum1M: =0; 
AbsErrSumEM: =0; 
j:=0; {initialize output 


END; {1F} 


PrintResults; 


{initialize absolute int errors} 


IF (fPrime=fMax-StepSize) THEN 


t NTE (fTwo):5:2,! a 
i‘ NIM(fIwol:5:2,’ 7, 
Nie 

7 INIEINIC oP ey) S's 2. ° ey. 


{integrate errors} 
AbsErrSumTE:=AbsErrSumTE+abs(ErrorTe); 


{initialize integrated errors} 


line counter} 


{int abs err} 
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BEGIN {IF} 
Writeln (' 4 


‘ ns 
: dc aca Ba etl 6) RS 
Writeln; 


Writeln (’ AVERAGE ERRORS: ‘, 
(Err sum) ey ji ico. 
LerrSumiIM/ 922722, 
GEnmSumE My Wlisede2 1 

Writeln (‘AVERAGE ABSOLUTE ERRORS: ‘, 
(AbsErrSumTE/j):12:2, 
(AbsErrSumTM/j):27:2, 
(AbsErrSumEM/j):27:2); 

ENDS sc Fal 


END; {OutputNvsf} 


PROCEDURE OutputDvsf (Method: 1..2); 


A A ts ee 


Tabulates the number of modes aluated spectral density 
the approximate formulae and from the differential form 
of the approximate formulae. 


VAR deltaTE,deltalTM,deltaEM: real; {values of D} 
PROCEDURE PrintTitles; 


iI 
{| Prints column headers for table. 
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BEGIN {PrintTitles} 


WritelIn (' at 
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Writeln PA Pet AT eit AU PREG 1 SUD 
END; {PrintTitles 


PROCEDURE PrintDifferences: 


{ 
{ 
{ 


Uses the difference relations to compute deltaN. 


BEGIN {PrintDifferences} 
deltaTE:=NTE(fTwo)-NTE(fOne) ; 
deltaTM:=NTM(fTwo)-NTM(fOne); 
deltaEM: =NEM( fTwo)-NEM(fOne) ; 
ErrorTE:=TE-deltalE; 
Error1TM:=TM-deltalM: 

ErrorEM: =EM-deltaEM; 


Wi tem (ale Omen) Oman 

reali (0 1 We etsy, | 7 olen eee, = 9 
prorvEesvsa,” 8 je 

Writes (7 \ariNeS.4 edelitalMes!2h ae .ae 
ERromiiNiaies 2. han) i 

Writes (“(77EM: 5,” ade)litac Mig 5 aun 
EGROReM is 72) ameein)o 

Writeln; 


ErrSumTE:=ErrSumTE+Errortve; 
AbsErrSumTE:=AbsErrSumTEt+abs(ErrorvE) ; 
ErrSumTM:=ErrSumTM+ErroriM:; 
AbsErrSumTM: =AbsErrSumTM+abs(Error TM) ;: 
ErrSumEM: =ErrSumEM+ErroreM; 
AbsErrSumEM: =AbsErrSumEM+abs(ErroreEM) ; 
hots Kase 
END; {PrintDif ferences} 


PROCEDURE PrintDifferentials;: 


{ 
{ 
{ 
{ 


Uses the differential relations to compute D. 


BEGIN {PrintDifferentials} 
deltaTE:=DTE(fOne,bandwidth) ; 
deltaTM:=DTM(fOne,bandwidth) ; 
deltaEM: =DEM( fOne,bandwidth) ; 
ErrorTE:=TE-deltaleE; 
Error™M:=1M-deltalIM; 

ErrorEM: =EM-deltaEM; 
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Writeln; 


ErrSumTE:=ErrSumTE+Errorle; 
AbsErrSum7TE:=AbsErrSumTE+abs(Errorve) ; 
ErrSum™M:=ErrSumTM+Error IM; 
AbsErrSum1TM: =AbsErrSumTMt+abs(Error TM) ; 
Err SumEM: =ErrSumEM+ErroreEM; 
AbsErr SumEM: =AbsErr SumEM+abs(ErroreEM) ; 
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K:=k+1; 


END; { 


BEGIN {Ou 
IF (#P 


BEG 


END 
IF (Me 
Reet 2 
BEG 


END; {Out 
BEGIN {Gener 


PrintDifferentials} 


tputDvsf } 
rime=fMin) THEN 


Tele} 
1F (Method=1) THEN 


BEGIN {IF} 
Writeln (' D versus Frequency Data ’, 
‘(Calculated from Differences):' ); 
Wri n; 
END { 


ELSE 


tel 
iE}, 


BEGIN {IF} 
Writeln (' D versus Frequency Data ', 
‘(Calculated from Differentials):'); 
Writeln; 
END saat 15} 


PrintTitles; 
ErrSumTE:=0; {initialize integrated errors} 
ErrSumTM:=0; 


ErrSumEM: =0; 

AbsErrSumTE:=0; {initialize absolute int errors} 
AbsErrSum™M:=0; 

AbsErrSumEM: =0; 

k:=0; {initialize output line counter} 


ee 


thod=1) THEN PrintDif ferences 
ELSE PrintDifferentials; 
rime=fMax-StepSize) THEN 


IN {15 ) 

Write ng ieee Saeed, 
bh tc ra ee eee 12) 

Writeln; 

Writeln ('’ AVERAGE ERRORS: ‘, 
(ErrSumTE/kK):21:2, 

(ErrSumTM/kK) :27:2, 
(ErrSumEM/kK) :27:2); 

Writeln (’ AVERAGE ABSOLUTE ERRORS: ‘, 
(AbsErrSumTE/kK):12:2, 
(AbsErrSumTM/k) :27:2, 
(AbsErrSumEM/kK) :27:2); 


putDvsf} 


ateData} 


IF (bandwidth=0) THEN 


BEGIN 
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fPrime:=fMin; {initialize fPrime} 
REPEAT 


fOne:=0; {set 1st marker} 
fTwo:=fPrime+StepSize; {set 2nd marker} 
TE:=0; {initialize TE mode counter} 

T™:=0; {initialize TM mode counter} 

EM:=0; {initialize total mode counter} 
CountModes(fOne,fTwo); {find N} 

OutputNvsf; {print values in a table} 
fPrime:=fPrime+StepSize; {increment fPrime} 


UNTIL fPrime>=fMax 


END {IF} 


BCoce OR 1. = lmOnee DO 


BEGIN {FOR} 


fPrime:=fMin;: {initialize fPrime} 
REPEAT 


fOne:=fPrime; {set 1st marker} 
fTwo:=fPrimet+bandwidth; {set 2nd marker } 
TE:=0; {initialize TE mode counter} 

T™:=0; {initialize TM mode counter} 

EM:=0; {initialize N} 
CountModes(fOne,fTwo); {find N} 
OutputDvsf(i); {print values in a table} 
fPrime:=fPrime+StepSize; {increment fPrime} 


UNTIL fPrime>=fMax 
END; 


{FOR} 


END; {GenerateData} 


BEGIN 


Readin (Lx,Ly,Lz, fMin, fMax, StepSize, bandwidth) ; 


wWriteln 


Writeln; 


Writeln 
writeln 
writeln 
Writeln 


Writeln:; 
ComoputeVSL; 


('The input parameters used in this run are:’:70); 


(exe yalzererG OP xc Oar ii ime Lyon men meme Zo meee ee 
(? aia tees SB) SCD AMO ia” tere, gtAiEOS Sg” ink2” I) 
(*>StepSize*s’ 742"Stepsize:5, "Hz 1)" 

(‘bandwidth =' :42,bandwidth:5,’ Hz’); 


GenerateData 
END. {ModeDistribution} 


SENTRY 
0.23 
0.50 
OF52 
4.0E8 

In SES 
5000000 
50000000 
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A.3 CARETAB 


10 

20 

30 

40 

50 

60 

70 

80 

90 

100 
110 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
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300 
B10 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
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450 
460 
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480 
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510 
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570 
580 
590 
600 
610 
620 


REM SSSSS SSS SS SSSSSSSseSS5SeSSeeSSSSSSSSeeeeeS5 5525555255252 5222522e2252==== 
ee F. PROGRAM: CARETAB VERSION: 1.1 REVISED: AUGUST 08, 1984 
REM SSSSS SS SSS SS SSS SSS SSeS SSSSSeeSSSee5eeSe5 SSeS S5 55S S25 eSe5 55225222 2ee=5=== 
REM PURPOSE: To calculate the conceivable number of TE, TM and EM 
REM modes in a (computed) series of rectangular cavities 

' REM within a given bandwidth. 
REM NOTES: The mode count results are sorted into four types: 
REM 1. TS, the total number of (mathematical) solutions 
REM to the dispersion relation. 
REM 2. NL, the total number of non-longitudinal solutions 
REM only (as given by MacKay, et.al.). 
REM 3. TE, the true number of Transverse Electric modes. 
REM 4. TE, the true number of Transverse Magnetic modes. 
REM Run with DOS Version 1.1 using the advanced interpreter. 
REM Allocate 32 Kbytes of RAM to the flashdisk. 
REM SSS SSS SS SS SSS SS SSS SSS SSS SSS SSS SS SSS SSS SSeS eS eee eeeeeeeeeeeeeeeSeeSE= 
REM SSSSSSS SSS SSS SSS SS SS SS SSS SS SSS SESS SESS eeEe Ee fo e2eeee2ee2eeeee2e2e2e2ee222e2Se=2=== 
ee Definitions and Dimensions (constants, variables, arrays) 


DEFSTR A ‘Character string (esp. user responses entered from keyboard) 
DEFSNG B,C,D,F,H,U,V,W,X,Y,Z ‘Single precision (real variables) 

DEFINT E,G,I,J,K,L,M,N,O,P,Q,R,S,T ‘integer (counters, integer variables) 
DIM CAVS(500,3) ‘contains dimensions of each cavity tested 

DIM MODS (500,5) ‘contains results of the various mode counts 

DIM MAXTE(100) ‘contains indices of CAVS and MODS (max. no. of TE modes) 
DIM MAXTM(100) ‘contains indices of CAVS and MODS (max. no. of TM modes) 
DIM MAXEM(100) ‘contains indices of CAVS and MODS (max. no. of EM modes) 
LET HZ=1000000! ‘MHz to Hz conversion factor 

LET C=3E+10 ‘speed of light in free space (cm/s) 

REM 
REM 
REM 


: Q=1 ‘indices used in arrays CAVS and MODS (for storing and printing) 
: T=0 ‘output device number and number of lines printed to that file 
0: RL=0 ‘logical variables used to direct program flow 

TEM=0: TMM=0: EMM=0 ‘reference mode counts used in MAXTE, MAXTM, MAXEM 
TEC=0: TMC=0: EMC=0 ‘counters and indices used in MAXTE, MAXTM, MAXEM 
OPEN "SCRN:" FOR OUTPUT AS #1 ‘establish CRT screen as output device 
OPEN "C:TEMP” FOR APPEND AS #2 ‘establish flash disk as output device 
OPEN "C:TEMP" FOR INPUT AS #3 ‘establish flash disk as input device 

OPEN "COM1:300,N,8" AS #4 ‘establish printer as output device via COM1 
ON KEY(1) GOSUB 1380 ‘interrupts program after completing current cavity 


KEY OFF: BEEP: CLS ‘erase soft key definitions and clear screen 

READEX1,X2 -XS7XINY 2, XS 7 olpee cov bia’ 2eteaGecefaultedata 

DATAS 200230 , 207207 L920) 2874490079807 "short-execution demo data 

DATA 48.5751. 59 0.5)38.5,4165,0.5,38.5,41.5,0.5,2425,2475,, sPuschner)s,data 
INPUT "New data? ", Al: IF NOT(LEFTS$(Al,1)="y") THEN 700 ‘or default data 
PRINT: PRINT TAB(20) "NOTE: Large bandwidths and/or ranges for X,Y,2" 
PRINT TAB(28) "will lead to VERY long execution times!": PRINT: PRINT 

REM some rudamentary error trapping is provided to catch likely typos 
INPUT "Enter FIRST target value for X/0 in cm...", Xl: IF X1[=0 THEN 550 
INPUT "Enter LAST target value for X|0 in cm...", X2: IF X2[=0 THEN 560 
IF (X1|X2) THEN PRINT: PRINT "X2 must not be less than Kose GOTO 50 
INPUT "Enter STEPSIZE (| 0) rope 6 TUN Nooo 7 OSG MISOUSAWR IG) NSS psoas, Syxeyo 
INPUT "Enter FIRST target value for Y|0 in cm...", Yl: IF Y1[=0 THEN 590 
INPUT "Enter LAST target value for Y|0O in Clinvecmeyey cite hey =O THENM OOO 
Ag! (Yas x2) THEN PRINT: PRINT "“Y2 must not be Jessmthaniny Juice GOLO S50 
INPUT "Enter STEPSIZE (| 0) LOL NE CMe mOsm PRUND el havo =O CHENG O20 
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CARETAB Cont'd, 


630 
640 
650 
660 
670 
680 
690 
700 
710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 
1110 
1120 
ls © 
1140 
1150 
1160 
T7.0 
1180 
1190 
1200 
1210 
1220 
1230 
1240 


INPUT "Enter FIRST target value for 2/0 in CMA Ga A AUG Aig Pal es) Uesbasl (ie 
INPUT "Enter LAST target value for 2/0 in com...", Z2: IF 22[=0 THEN 640 
IF (21122) THEN PRINT: PRINT "Z2 must not be less than 21!": GOTO 550 
INPUT "Enter STEPSIZE (|0) for Z in cm...", ZS: PRINT: IF ZS[=0 THEN 660 
INPUT "Enter LOWEST frequency (=|0) of BW in tion crs Jwala igo) Igbl io) Beebo) (3 7/(0) 
INPUT "Enter HIGHEST frequency (|0) of BW in MHz...", F2 

IF (Fl|F2) THEN PRINT: PRINT "Fl must be greater than F2!": GOTO 670 

CLS: LOCATE 25,10: COLOR 0,7 ‘position cursor and change to reverse video 
PRINT " Press Fl to interrupt program after computing next cavity. "; 
COLOR 7,0: PRINT: LOCATE 1,1: BEEP ‘reset font and reposition cursor 


FOR X=Xl1 TO X2 STEP XS: FOR Y=Yl TO Y2 STEP YS 'X and Y dimension loops 
REM space 1] line and print column headers every 2S rows 

PRINT: PRINT TAB(2) “X(cm)" TAB(10) "Y(cm)" TAB(18) "Z(cm)"; 

PRINT TAB(31) "TS" TAB(37) "NL" TAB(43) "TE" TAB(49) "TM" TAB(55) "EM" 
PRINT: PRINT USING "###.## "; X;Y; ‘print the X and Y dimensions 

FOR Z=Z1 TO Z2 STEP ZS 'Z dimension loop 

CAVS (P,1)=X: CAVS(P,2)=Y: CAVS(P,3)=Z 'store the current cavity dimensions 
PRINT TAB(17);: PRINT USING "###. #44 "; Z; ‘print the 2 dimension 
KEY(1) ON: KEY(1) STOP 'trap Fl and execute before starting next loop 


TS=0: NL=0: TE=0: TM=0: TIMES="00" ‘reset all mode counters and timer 
FOR L=0 TO FIX (X*2*F2*HZ/C)+1 STEP 1 ‘first mode index 

FOR M=0 TO FIX (Y*2*F2*HZ/C)+1 STEP 1 ‘second mode index 

FOR N=0 TO FIX(2Z*2*F2*HZ/C)+1 STEP 1 ‘third mode index 

F= (C/2) *SQR(((L/X)~2)+ ( (M/Y)~2)+((N/Z)~2))/HZ ‘resonant frequency in MHz 


REM count the total number of (mathematical) solutions to equation for F 
IF (F([Fl) OR (F|F2) THEN 1030 ELSE TS=TS+1 ‘ignore if F is out of range 
REM count the non-longitudinal modes (as Mackay, et al.) 
IF (L{O AND N|O) OR (M|O AND N|0) OR (L|O AND M|0) THEN NL=NL+1 
REM count TE and TM modes separately and combine them outside the loop 
IF (L|O AND N|O) OR (M|O AND N|0) THEN TE=TE+1 ‘count the TE modes 
IF (L|O AND M|0) THEN TM=TM+1 ‘count the TM modes 
NEXT N,M,L 
EM=TE+TM 'comptte the total number of electromagnetic modes 


IF TE[TEM THEN 1150 ‘continue if less than current maximum 

IF TE=TEM THEN 1100 ELSE TEC=0 ‘test for a tie and reset counter if not 
TEC=TEC+1: TEM=TE ‘increment counter and establish TE as new maximum 
MAXTE (TEC)=P ‘enter corresponding CAVS and MODS index in array MAXTE 


IF T[TMM THEN 1220 ‘continue if less than current maximum 

IF TM=TMM THEN 1170 ELSE TMC=0 ‘test for a tie and reset counter if not 
TC=TMC+1: TMM=TM ‘increment counter and establish TM as new maximum 
MAXTM (TMC)=P ‘enter corresponding CAVS and MODS index in array MAXTM 


REM LEE bad deech: dated aight eee Sak ad sean + peer et SER Coot ae Cee a Sra === === 
IF (TE+TM) [EMM THEN 1290 ‘continue if less than current maximum 


IF EM=EMM THEN 1240 ELSE EMC=0 'test for a tie and reset counter if not 
EMC=EMC+1: EMM=EM ‘increment counter and establish EM as new maximum 
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CARETAB Cont'd. 


1250 
1260 
1270 
1280 
1290 
1300 
L320 
1320 
1330 
1340 
1350 
1360 
1370 
1380 
1390 
1400 
1410 
1420 
1430 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
5310 
1540 
1550 
1560 
1570 
1580 
590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
TO 
1780 
1790 
1800 
1810 
1820 
1830 
1840 
1850 
1860 


eee nd "enter corresponding CAVS and MODS index in array MAXEM 


REM Storing and Printing of Results to the Screen 

MODS (P,1)=TS: MODS (P,2)=NL: MODS (P,3)=TE: MODS (P,4)=TM: MODS (P,5)=EM 

PRINT TAB(29);: PRINT USING "#### "; TS,NL,TE,TM,EM; ‘print mode counts 
IF Z=21 THEN PRINT "* ";TIMES ELSE PRINT 'print execution time once 
P=P+1 ‘increment array index (cavity counter) 

NE XD Zig y 9X 

REM SS SSS SS SSS SS SS SS SSS SS SSS SSS SSS SSS SSS SSS SS S222 e225 eS e2e 252 2Ee2e2e5EEe5SEeeEee===== 
REM (Optional) Storage of Mode Count Results in file C:TEMP 

REM disable trapping and ask whether or not storage is required 


GOSUB 1960: IF PL THEN GOTO 1540 ELSE S=2 ‘device #2 is file TEMP 

FOR I=1 TO P-1 STEP 1+FIX((Z2-21)/ZS) ‘intervals for column headers 
PRINT#S, ‘blank line ; 

PRINT#S,TAB(12) "X(cm)" TAB(20) "Y(cm)" TAB(28) "Z(cm)"; ‘column headers 


PRINT#S,TAB(41) "TS" TAB(47) "NL" TAB(53) "TE" TAB(59) "TM" TAB(65) "EM" 
PRINT#S,: PRINT#S,TAB(11); ‘blank line, position cursor 

PRINT#S,USING "###.## "; CAVS(Q,1);CAVS(Q,2); 'print X and Y 

FOR J=0 TO FIX((2Z2-21)/ZS) STEP 1 ‘line items per interval 
PRINT#S,TAB(27): PRINT#S,USING “###. ## "5 GNS(O_3)) 7 Vjskeaine ve 
PRINT#S,TAB(39); ‘tab and print mode count results 

PRINT#S,USING “#### "; MODS(Q,1),MODS(Q,2) ,MODS(Q,3) ,MODS(Q,4) ,MODS (Q,5) 


Q=Q+1 ‘increment array index 

NEXT J,1I 

REM SS 4 
REM Printing of Cavities Identified Above 

REM — ee ee 
S=l1 'reset output device switch S to screen 


PRINT#S,"The cavities which support the most TE modes are:": GOSUB 2310 


PRINT#S,: GOSUB 2310 ‘call PAGEPAUSE to check line count after each line 
FOR I=1 TO TEC STEP 1 ‘print cavities with the most TE modes 

PRINT#S, TAB(5); ‘tab and print dimensions of TE cavities 

PRINT#S," X=";: PRINT#S,USING "###.##"; CAVS (MAXTE(I),1);: PRINT#S," cm"; 
PRINT#S,” Y=";: PRINT#S,USING "###.##"; CAVS (MAXTE(1I),2);: PRINT#S," cm"; 
PRINT#S," Z=";: PRINT#S,USING "###.##"; CAVS (MAXTE(1),3);: PRINT#S," cm" 
GOSUB 2310: NEXT I 

PRINT#S,: GOSUB 2310: PRINT#S," Each cavity supports"; 

PRINT#S,USING "####"; MODS (MAXTE (TEC) ,3);: PRINT#S," TE modes." 


GOSUB 2310: PRINT#S,: GOSUB 2310: PRINT#S,: GOSUB 2310 'blank lines 
PRINT#S,"The cavities which support the most TM modes are:": GOSUB 2310 
PRINT#sS,: GOSUB 2310 


FOR I=l TO TMC STEP 1 ‘print cavities with the most TM modes 

PRINT#S, TAB(5); ‘tab and print dimensions of T™ cavities 

PRINT#S,"” X=";: PRINT#S,USING "###.##"; CAVS (MAXTM(I),1);: PRON T4 Siac iniaas 
PRINT#S," Y=";: PRINT#S,USING Pees CAV SI(MAXTIM (le)i 2); tee RUNG tone CIN arr 
PRINT#S," Z=";: PRINT#S,USING "##2.9#">; CAVS (MAXTM(I),3);: PRINT#S," cm" 
GOSUB 2310: NEXT I 

PRINT#S,: GOSUB 2310: PRINT#S,” Each cavity supports"; 

PRINT#S,USING "####"; MODS (MAXTM (TMC) ,4);: PRINT#S," TM modes." 

GOSUB 2310: PRINT#S,: GOSUB 2310: PRINT#S,: GOSUB 2310 ‘blank lines 


PRINT#S,"The cavities which support the most EM modes are:": GOSUB 2310 
PRINT#S,: GOSUB 2310 

FOR I=1 TO EMC STEP 1 'print cavities with the most EM modes 

PRINT#S, TAB(5); ‘tab and print dimensions of EM cavities 


PRINT#S," X=";: PRINT#S,USING "###.##"; CAVS (MAXEM(I),1);: PRIND$ Sic Min 
PRINT#S," Y=";: PRINT#S,USING "38# #8": CAVS (MAXEM(I),2);: PRINT#S," cm"; 
PRINT#S," Z="73: PRINT#S,USING "###.##"; CAVS (MAXEM (1) ,3)7: PRINT#S," cm" 
GOSUB 2310: NEXT I i 
PRINT#S,: GOSUB 2310: PRINT#S," Each cavity supports”; 


PRINT#S,USING "####"; MODS (MAXEM (EMC) ,5);: PRINT#S," EM modes." 
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CARETAB Cont'd. 


1870 
1880 
1890 
1900 
1910 
1920 
1930 
1940 
1950 
1960 
OHO 
1980 
1990 
2000 
2010 
2020 
2030 
2040 
2050 
2060 
2070 
2080 
2090 
2100 
2110 
2120 
213 0 
2140 
2150 
2160 
2170 
2180 
2190 
2200 
2210 
2220 
2230 
2240 
2250 
2260 
2270 
2280 
2290 
2300 
2310 
2320 
2330 
2340 
2350 
2360 
2370 
2380 
2390 
2400 
2410 
2420 


GOSUB 2310: PRINT#S,: GOSUB 2310: PRINT#S,: GOSUB 2310 ‘blank lines 

IF S=2 THEN 1900 ELSE GOSUB 1970 ‘ask if storage in file TEMP is required 
IF PL THEN 1900 ELSE S=2: GOTO 1550 ‘rewrite results to file C:TEMP 

GOSUB 2100 

nie ON ‘replace soft key definitions 


KEY(1) OFF ‘disable Fl trapping on first call 

BEEP: LOCATE 25,1: PRINT STRINGS (80,0); 'get attention, clear prompt line 
LOCATE 25,14: COLOR 0,7 'position cursor and change to reverse video font 
PRINT " Do you want these results saved in the file C:TEMP? "; ‘prompt 
COLOR 7,0 ‘return to normal video font 

INPUT " ", A2: A2=LEFTS$(A2,1) ‘examine first character of response A2 
LOCATE 25,1: PRINT STRINGS (80,0); ‘clear prompt line 

IF (A2="Y") OR (A2="y") THEN PL=0: RL=1: RETURN 'storage required 

IF (A2="N") OR (A2="n") THEN PL=l]: RETURN ‘storage not required 

PRINT "Please answer YES or NO (Y or N).": PRINT: GOTO 1970 ‘instructions 
RETURN 


IF RL=0 THEN RETURN 'check whether or not any results have been stored 
BEEP: LOCATE 25,1: PRINT STRINGS (80,0); ‘get attention, clear prompt line 
LOCATE 25,14: COLOR 0,7 "position cursor and change to reverse video font 
PRINT " Do you want these results sent to the printer? "; 

COLOR 7,0 ‘return to normal video font 

INPUT " ", A2: A2=LEFTS (A2,1) ‘examine first character of response A2 
LOCATE 25,1: PRINT STRINGS (80,0); ‘clear prompt line 

IF (A2="Y") OR (A2="y") THEN S=4: GOTO 2200 ‘set device switch to printer 
IF (A2="N") OR (A2="n") THEN RETURN 'continue if no printing required 
PRINT "Please answer YES or NO (Y or N).": PRINT: GOTO 2110 ‘instructions 
BEEP: LOCATE 25,1: PRINT STRINGS (80,0); ‘get attention, clear prompt line 
LOCATE 25,14: COLOR 0,7 ‘position cursor and change to reverse video font 
PRINT " Turn the printer on and strike the SPACE BAR to Start printing. "; 
A3=INKEYS: IF NOT (A3=CHRS (32)) THEN 2230 'pause until space bar is struck 
COLOR 7,0: PRINT; ‘return to normal video font 

LOCATE 25,1: PRINT STRINGS (80,0); ‘clear prompt line 

WHILE NOT EOF(3): LINE INPUT#3,A4: PRINT#S,A4: GOSUB 2310: WEND 

RETURN 


REM ee 
IF S=l AND T=|21 THEN 2340 'pause when CRT screen is full 

IF S=2 OR S=3 THEN RETURN ‘ignore call 

IF S=4 AND T=|53 THEN PRINT#S,STRINGS (11,10): GOTO 2340 ELSE T=T+1: RETURN 
BEEP: LOCATE 25,1: PRINT STRINGS (80,0); ‘beep attention, clear prompt line 
LOCATE 25,22: COLOR 0,7 ‘position cursor and change to reverse video font 
PRINT " Strike the SPACE BAR to continue. "; ‘prompt 

A5=INKEYS: IF NOT(A5=CHR$(32)) THEN 2370 'pause until space bar is struck 
COLOR 7,0: PRINT ‘return to normal video font 

LOCATE 25,1: PRINT STRINGS (80,0) ‘clear prompt line 

T=0 'reset line counter 

RETURN 

END 
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APPENDIX B. CAVITY MODE COUNT REFERENCE TABLES 


Table B.1 Numerically determined mode count reference data 
for several rectangular cavities excited by frequencies 
within the range 2425 < £ < 2475 MHz. TS is the total 
number of mathematical solutions to the dispersion relation; 
NL is the number of nonlongitudinal modes; TE and T are the 
actual numbers of TE and TM modes respectively; EM is the 


total number of electromagnetic modes. 
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Table B.2 Numerically determined mode count reference data 
for several rectangular cavities excited by frequencies 
within the range 900 < f < 930 MHz. TS is the total number 
of mathematical solutions to the dispersion relation: NL is 
the number of nonlongitudinal modes; TE and TM are the 
actual numbers of TE and TM modes respectively; EM is the 


total number of electromagnetic modes. 
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Table B.3 A comparison of 


several cavities of typical sizes used 


ovens. 
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Table B.4 Variations in the mode count data presented in 


Tables B.1,2,3. 


BANDWI DTH RANGE OF CAVITY DIMENSIONS RANGE OF MODE COUNTS 
(MHz) X (cm) Y (cm) Z (cm) Ts NL TE ™ EM 
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APPENDIX C. ON THE CROSS-COUPLING PROBLEM 


Solid state microwave power applicators have been 
technologically feasible for over a decade. However, in 
order to make them economically feasible (given the present 
or slightly advanced state of the art), efficient means must 
be devised for combining the outputs of several medium power 
microwave sources. The simplest method is to use the cavity 
itself as the power combining network, by separately 
coupling each source to the cavity through an independent 
feed Structure. The Mcross-coupling® probleme Mi s@ethus 


defined as achieving this independence. 


Several methods have been hypothesized which should, 
theoretically, reduce the cross-coupling between multiple 
sources. For example, operating each source at widely 
different frequencies (e.g. 915 vs. 2450 MHz) or using 
cross-polarized antennae are commonly cited [48]. However, 
to the author's knowledge, no one has examined the problem 
experimentally to determine the severity (or nonseverity) of 
thissproblem qualitatively.™ Part off—the reason¥efor @this 
lack of information may be the difficulty in generalizing 
any experimental data. And for this reason, the remainder 
of the discussion which follows is primarily conjecture, but 


is based on observations made while conducting’ the 


measurements of Chapter 4. 
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The measurement system described in Chapter 4 is easily 
modified to observe the cross-coupling between any two 
antennae. Basically, a transmission loss measurement can be 
performed by connecting the output of a receiving antenna to 
channel B of the swept amplitude analyzer through a _ second 


detector as indicated in Figure 4.1. 


It was found that when the receiving probe was very 
short, and the cavity was empty, only cross-coupling through 
normal cavity eigenmodes can occur. In this case, the 
transmission spectrum is discrete and differs mainly in 
magnitude from the input spectrum. The degree of 
cross-coupling for some of the modes within the domain of 
Figure 4.2 were classified as strong (S), medium (M) or weak 
(WieeraSste indicated —an (Figure Cale The actual level of 


transmission loss was small. 


When the receiving probe waS a realistic length 
(~10 cm), and a mode stirrer was included, the transmission 
spectrum was observed to be virtually continuous, and the 
transmission loss was much higher than in the above case. 
The continuous spectrum can probably be explained as a 
result of the length and proximity of the two antennae which 
allowed cross-coupling through near fields or evanescent 
modes (which is enhanced by the mode stirrer). Theres 
observation would seem to indicate that the cross-coupling 


problem is indeed severe ina practical case. However, many 
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Figure C.1 Measured spectral distribution of modes (NO) 
for the experimental cavity discussed in Chapters 4 and 5, 
showing the relative cross-coupling between top and side 


mounted antennae. The receiving antenna is less than 1 cm 


in length. 
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new power transistors can withstand infinite VSWR in a class 
C mode (at least temporarily). And when the cavity is 
heavily loaded with lossy dielectric, the actual level of 
received power due to the cross-coupling effect would 
presumably be much less than the output capability of the 
source. In any case, the output spectrum of solid state 
power amplifiers is broad enough that much _ power is 
neflectedé,from Eangremptyancavity sevens whene thefinput is 
perfectly matched to the cavity (at the eigenfrequencies). 
Hence, although the cross-coupling per se can be serious, it 
may not necesSarily be detrimental to the operation of a 
solid state applicator. More thorough experimental studies 
with transistor sources are required in order to determine 


the true extent of this problem. 


Numerical Models 


Since general results are difficult to obtain, either 
experimentally or analytically, it would be extremely useful 
et a numerical algorithm were available that could 
parametrically model the cross-coupling between two (or 
more) antennae. Again, a completely general model could not 
be obtained with a finite number of parameters, but perhaps 


certain classes of antennae could be modeled. 


It may be possible to base such a model on a 


perturbation-type analysis, where the perturbed 
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eigenfunctions are numerically integrated over the receiving 
antenna. However, some philosophical problems regarding the 
applicability of perturbation theory to a cavity which is 
not a closed thermodynamic system, may exist. A more 
complicated approach, but one which is (probably) free of 
philosophical questions, is to employ the image theory of 
Kelvin. That is, once all of the image antennae are 
determined, the theory of antennae arrays could perhaps be 
used to compute the image currents on the receiving antennae 
due to the exciting antennae. If this method is feasible, 
it may be the most general since free space radiation 
patterns for arbitrarily shaped antennae can be determined 


empirically. 
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APPENDIX D. ASYMPTOTIC FORMULAE FOR WAVEGUIDES 


Recall that the two dimensional EM boundary value 
problem leads to an eigenvalue equation for the cutoff 
frequencies of waveguide modes. Asymptotic formulae for the 
spectral distribution and density of these eigenfrequencies 
can be derived in exactly the same manner as that used in 
Chapter 5 to obtain the asymptotic formulae for cavity 


modes. 


The number of characteristic points in a two 
dimensional p-space, to first order, is just the area 
enclosed in the first quadrant by the ellipse described by 


EGe (2.29 )EViz. 


=A es ee NO (w) (Dip 
Arm Ce 


where A is the cross-sectional area of the guide (normal to 
the x3; axis). N¥(w) is thus the two dimensional counterpart 
of No(w), the ambivalent first term of the asymptotic 


expansion for cavity modes. 


Correcting ssLor the undercounting/overcounting of 


characteristic points on the coordinate axes of p-space, one 


obtains 


Nrylw) = No(w) - (Li+b,) wo + 1 (Diz) 
27 c g 
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for the TM modes, and 


Nye(w) = NO(w) + (Ly+L2) w - 3 (Desh 
27 Cc 4 


for the TE modes. Thus the total number of EM modes with 


cutoff frequencies not greater than w is approximately 


Nen(w) = 2N¥ (w) oad a (D4) 
2 


and the total number of EM modes with cutoff frequencies 


between w and wtdw 1s approximately 


Dew) dw = A dw = 2D (O) 60 CDie5y) 
T 


WwW 
Cc? 


where D¥(w) is the two dimensional counterpart to Do(w) as 
defined in Chapter 5. Hence, the spectral density function 


for waveguide modes is 
Dem(w) = 2D¢(w) (D.6) 


which is just the derivative of Ne,(w) with respect to w. 
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